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1.  INTRODUCTION 

In  recent  years  the  Marine  Physical  Laboratory  of 
Scripps  Institution  of  Oceanography  has  developed  and 
fielded  several  imaging  systems  for  assessment  of  the 
cloud  field  dynamics,  as  well  as  for  visibility  determina¬ 
tion.  The  Whole  Sky  Imager  (WSI)  images  the  sky  dome 
and  identifies  cloud  distributions,  from  which  earth-to- 
space  cloud  free  line  of  sight,  cloud  free  arc,  and  cloud 
cover  determinations  may  be  made.  The  Horizon  Scan¬ 
ning  Imager  (HSI)  images  scenes  near  the  horizon,  with 
a  small  field  of  view  lens  and  near-photopic  filtering,  for 
determining  the  visibility  over  extended  paths. 

Both  of  these  units  have  been  combined  and  adapted 
for  near-real  time  assessments,  in  the  Automated  Ob¬ 
serving  System  (AOS)  stationed  at  Otis  ANGB.  This 
woiic  has  been  reported  previously  in  GL-TR-89-0061, 
Johnson,  et  al,  (1989),  and  PL-TR-9 1-201 6(1),  Johnson, 
et  al  (1990),  and  PL-TR-9 1-22 16,  Johnson,  et  al  (1991). 

During  the  past  year,  work  has  continued  to  upgrade 
the  performance  of  the  daytime  HSI  system,  and  extend 
the  HSI  capability  to  the  nighttime  regime.  This  final 
report  discusses  this  work,  as  well  as  summarizing  rec¬ 
ommended  upgrades  to  the  AOS  system. 

2.  THE  COMPOSITE  HSI/WSI  OBSERVING 
STATION 

The  composite  HSI/WSI  observing  station  includes  a 
daytime-capable  HSI,  and  day-time  capable  WSI.  In 
addition,  the  HSI  hardware  and  software  capabilities 
have  been  extended  to  be  night-capable,  i.e.  provide 
capability  for  determining  visibility  at  night,  although 
the  night  capability  is  not  fully  automated  at  this  time.  A 
night-capable  WSI  unit,  for  measurement  of  cloud  fields 
at  night,  has  been  developed  under  separate  contract,  and 
is  a  reasonable  candidate  for  inclusion  in  the  observing 
station  at  a  later  point  in  time. 

This  section  will  provide  a  summary  overview  of  the 
daytime  HSI  and  WSI  components.  The  Night  HSI  will 
be  discussed  further  in  later  sections. 

2.1  The  Daytime  HSI 

The  Horizon  Scanning  Imager  uses  a  sol  id  state  Cha  rge 
Injection  Device  (CID  sensor  for  the  acquisition  of  5 1 2 
X  480  images  of  the  horizon.  The  use  of  a  5  degree  field 
of  view  lens,  in  combination  with  an  auto-iris  and  a 
spectral  filter  (green,  with  a  peak  wavelength  of  550  nm, 
and  a  passband  of  70  nm),  allows  the  acquisition  of  high 
resolution  near-photopic  images.  A  precision  rotary 
table  is  used  for  highly  reproducible  scene  acquisition 
around  the  horizon. 


Dark  targets  arc  previously  identified  in  the  scene  at  a 
variety  of  rangc.s.  A  .sampte  visibility  image  is  shown  in 
Fig.  1 .  In  this  scene,  the  relative  radiance  of  each  of  these 
targets  is  determined,  along  with  the  relative  radiance  of 
the  horizon  at  tlic  same  azimuth.  From  these  dclcmiina- 
tions,  the  apparent  contra.st  of  each  target  with  respect  to 
the  horizon,  Q,  is  determined. 

C,  =  (2.1) 

bl-r 

where 

tLr  =  Apparent  target  radiance 
bLr  =  Apparent  background  radiance 


Fig.  1  Sample  Visibility  Image.  Image  was  acquired  with 
the  HSI,  and  resulting  visibility  computations  are 
performed  automatically  in  the  field. 

Figure  2  illustrates  the  decrease  in  contrast  of  a  black 
target  with  respect  to  the  horizon,  as  range  to  the  target 
increases.  Note  that  the  apparent  contrast,  i.e.  the  con¬ 
trast  as  measured  from  a  given  range,  decreases  expo¬ 
nentially.  This  rate  of  decrease  depends  on  the  attenua¬ 
tion  coefficient  of  the  atmosphere.  Visibility  may  be 
defined  as  the  range  at  which  this  contrast  decreases  to 
the  human  contrast  threshold,  taken  to  be  .05  in  this 
illu.stration. 

The  lowercurvcs  of  Fig.  2  shows  the  apparent  contrast 
of  sample  targets  which  are  not  perfectly  black  (inherent 
radiance  not  equal  to  - 1 ).  The  HSI  is  able  to  measure  the 
contrast  of  non-ideal  tarets  such  as  represented  by  the 
lower  curves,  at  a  variety  of  ranges,  and  mathematically 
determine  the  visibility.  The  resulting  equation  dis¬ 
cussed  in  detail  in  Johnson,  et  al,  (1989),  is  given  below. 


V  =  r 


(2.2) 


This  technique  has  several  advantages  over  visual 
determination  of  the  visibility.  It  can  use  targets  which 
are  not  black,  which  are  not  at  a  range  equal  to  visibility, 
and  which  are  not  physically  adjacent  to  the  horizon.  In 
common  practice,  the  weather  observer  is  forced  to  use 
targets  with  these  limitations,  to  varying  degrees,  but  is 
unable  to  rigorously  compensate  for  the  non-ideal  condi¬ 
tions. 


Figure  2  Apparent  Contrast  vs.  Range.  Visibility  may  be 

defined  as  the  range  at  which  the  apparent  contrast 
of  the  ideal  target  drops  to  the  human  visual 
threshold,  shown  as  .05  in  this  illustration. 


In  comparison  with  point  scattcrmctcrs  andtransmis- 
someters,  the  HSI  has  the  advantage  of  using  a  long  path 
to  the  target,  which  yields  a  much  more  representative 
determination.  The  ability  to  measure  sector  visibility, 
i.e.  the  variance  around  the  horizon,  is  also  important  for 
many  applications. 

The  basic  functional  characteristics  of  the  HSI  system 
have  remained  essentially  as  described  in  John.son,  et  al, 
(1989).  The  general  concept  of  its  operational  sequence 
is  illustrated  in  Fig.  3.  The  diagram  has  of  course  been 
simplified  for  clarity  and  docs  not  include  the  myriad  of 
housekeeping  and  control  details  that  enable  the  auto¬ 
matic  image  acquisition  .sequences. 

The  HSI  continues  to  operate  reliably,  with  reason¬ 
able  accuracy.  This  year’s  upgrades  to  the  daytime 
visibility  capability  arc  di.scusscd  in  Section  4. 

2.2  The  Daytime  WSI 

Like  the  HSI,  the  Whole  Sky  Imager  is  a  passive, 
ground-based  sensor.  It  uses  a  fisheye  lens  to  acquire 
images  of  the  upper  hemisphere  down  to  an  80  degree 
zenith  angle.  A  combination  of  spectral  and  neutral 
density  filters,  in  conjunction  with  a  (5lD  camera,  enables 
acquisition  of  images  which  may  be  fully  calibrated 
radiometrically.  The  512  x  480  image  resolution  yields 
a  1/3  degree  spatial  resolution  in  object  space,  for  a 
footprint  of  roughly  1 7  meters  for  a  cloud  layer  at  3  km 
altitude. 


Fig,  3  HSI  Procedural  Flow  Chart 


The  cloud  data  undergo  several  processing  steps  to 
yield  a  final  cloud  decision  image,  in  whicli  areas  identi¬ 
fied  as  opaque  cloud  arc  displayed  as  white,  thin  cloud  arc 
displayed  as  yellow,  and  no  cloud  (clear)  arc  displayed 
with  blue.  Data  whichare  not  validduclocondilions  such 
as  offscale  dark,  are  also  identified  and  so  labeled.  Ir  this 
procedure,  a  number  of  calibrations  are  applied,  to  ac¬ 
count  for  slight  sensor  non-linearities,  filter  passband 
differences,  neutral  density  spectral  bias,  and  so  on. 
Images  acquired  at  650  nm  (red)  and  450  nm  (blue)  are 
ratioed  to  remove  most  of  the  directional  and  temporal 
variation  in  the  background  sky  radiance.  In  the  remain¬ 
ing  ratio  image,  opaque  and  thin  clouds  arc  identified  by 
ratio  thresholds,  yielding  a  determination  such  as  shown 
in  Fig.  4. 

This  ratio  technique  more  fully  illustrated  in  Fig.  5, 
allows  a  fully  automated  determination  at  each  pixel 
location  in  the  scene.  Itissignificantly  more  reliable  than 
identification  based  on  the  radiance  scene,  both  due  to  tltc 
removal  of  the  directional  bias  mentioned  above,  and  also 
because  clouds  can  often  be  darker  than  adjacent  sky 
regions. 

A  further  sophistication  to  the  cloud  identification 
scheme  uses  an  angularly  dependent  clear  sky  back¬ 
ground  ratio,  normalized  by  aerosol  load  or  haziness,  to 
correct  forthe  residual  directional  bias  in  the  ratio  image. 
This  technique,  developed  under  separate  contract,  has 
not  been  applied  to  the  composite  HSfAVSI  system.  Its 
advantages  will  be  discussed  in  Section  5. 


rijj.  4  Cluiui  DL’ci.sinn  Image.  Each  pixel  Itvaiion  ha.x  been 

uienliiiei)  iis  opaijiie,  ihin,  or  no  eloiul  by  ;in  .luloinaieti 
algiuithni  applied  lo  WSl  imagery. 


2.3  The  As-Built  Composite  HSI/WSI 

The  a.s-buill  HSl/WSI  unit  is  illustrated  in  Fig.  6, 
'lltese  tw(^  separate  cameia  units  are  controlled  by  an 
IBM  AT-cionemicrocomiJU  ter  with  an  image  board,  and 
a  2.2  Gbyte  tape  system  for  data  archival.  An  accessory 
control  panel,  illustrated  in  Fig7,allowscilhermanuaIor 
computer  control  of  the  various  functions  such  as  filter 


■ASlSIMAOenr  CORHEtTEPIWAfim  E2ttES2SiItBAIl,5  EEUyiflASLtEATAEAfE 


Fig.  5  WSI  Basic  Image  PrtKc.ssing  lioxv  C’fuiri 
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changer  selection,  solar  occultor  angle,  rotary  table  azi¬ 
muth,  and  so  on. 


Fig.  6  Composite  HSIAVSI  As-built  Field  Installation 


Fig.  7  Composite  HSIAVSI  Accessory  Control  Panel 


the  composite  HSIAVSI  system  the  two  sub-assemblies 
may  be  operated  jointly  under  an  operator-specified 
time-share  duty  cycle. 

A  nominal  automatic  duty  cycle  is  for  the  HSI  to  run 
horizon  sweeps  each  minute  for  eight  minutes  and  every 
ten  minutes  in.sert  a  cloud  cover  measurement.  In  this 
mode,  the  system  will  display  its  output  products  of 
sector  visibilities  and  cloud  cover  to  the  CRT  display 
and/or  output  the  data  to  printer  forhardcopy.  It  will  also, 
at  the  operator’s  discretion,  archive  both  original  imag¬ 
ery  and  derived  numerical  products  to  its  Exabyte  tape 
sub-system  for  later  retrieval  and  analysis. 

3.  VISIBILITY  AT  NIGHT 

This  .section  discusses  the  extension  of  the  H.SI  capa¬ 
bilities  to  enable  visibility  determinations  at  night.  First 
the  theory  of  visibility  at  night,  and  its  implication  for  the 
HSI  algorithms,  are  discussed.  The  revision  of  the  HS! 
hardware  is  then  described.  The  capabilities  and  limita¬ 
tions  of  the  CID  camera  for  night  time  use  arc  discussed, 
followed  by  sample  calculations. 

3.1  Definitions  of  Visibility  at  Night 

Douglas  and  Booker  (1977)  provide  an  in-depth  dis¬ 
cussion  of  visibility  at  night;  a  more  succinct  review  is 
given  in  Gordon  (1989).  This  section  gives  a  summary 
of  the  results. 

3.1.1  Basis  for  the  defining  equation  used  with  tran.s- 
missometers 

The  visibility  at  night  has  historically  been  defined  in 
terms  of  point  sources.  A  point  source  with  intensity  I,  in 
a  perfectly  transmitting  atmosphere,  will  produce  an 
illuminance  E  given  by 

E  =  l/r2  (3.1) 

where  I  is  given  in  candles,  and  E  in  lumen/km^,  or  km- 
candle  (Douglas  &  Booker,  p.  2-5).  Note  that  even  in 
clear  air,  E  is  a  function  of  range.  (The  luminance  or 
radiance  of  extended  sources  do  not  decrease  as  a  func¬ 
tion  of  r,  in  the  absence  of  attenuation.) 

In  a  real  atmosphere,  there  is  attenuation  which  can  be 
characterized  by  the  beam  transmittance,  defined  as  the 
ratio  of  transmitted  to  incident  flux.  The  illuminance  at 
range  r  becomes 


The  prototype  visibility  and  cloud  detection  systems 
were  fabricated  as  two  separate  hardwarc/softwarc  enti¬ 
ties  which  need  not  necessarily  be  integrated  into  a  single 
electro-mechanical  unit.  Each  may  stand  alone  and 
function  independently,  or  i  f  desired ,  under  the  control  of 


E  =  T,  I/r^  (3.2) 

were  T^  is  uniilcss.  The  attenuation  of  the  atmosphere 
can  be  equally  characterized  by  the  transmissivity  T 
(which  is  the  transmittance  over  unit  distance),  or  the 
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extinction  coclTicicnt  o  (which  is  the  nilo  ol  chat)gc  of 
the  flux).  These  icmis  arc  related  to  transmittance  by 

T,  :  r  -  c  .  (.T^i 

The  Federal  Meteorological  Handbook  (1988)  de¬ 
fines  tlic  visibility  as  “the  greatest  distance  at  which 
selected  objects  can  be  seen  and  idetiiiricd,"  whet,  the 
visibilit)'  markers  arc  defined  as  “dark  or  nearly  dark 
objects  view  cd  against  the  hori/on  sky  during  die  day,  or 
unfocused  lights  of  moderate  intensity  (about  25  cd) 
during  the  night.” 

This  is  a  psychophysical  definition.  Thus  at  night, 
using  Eq.  3.3,  the  definition  of  night  visibility  might  be 
stated  as  the  distance  at  which  the  illuminance  reaches  a 
certain  psychophysically  defined  llireshold.  when  the 
source  intensity  is  25  candles. 

Let  night  visibility  V  be  the  distance  at  which  the 
illuminance  reaches  a  threshold  E(.  Then  letlii  g  the 
range  r  =  V  in  Eq.  3.2,  and  using  Eq.  3.3.  the  defining 
equation  for  night  visibility  becomes 

E,  =  I7''/V-  (3.4) 

As  discussed  in  Douglas  and  Booker  P.  4-16,  this  is  the 
most  common  form  of  Allard’s  law,  and  may  be  consid¬ 
ered  as  the  defining  equation  for  vi.sibilit y  at  night. 

The  problem  with  this  definition  is  that  Et  itself  is  not 
well  behaved.  There  is  considerable  discussion  in  Dou¬ 
glas  and  Bookcrof  die  best  value  to  u.se,  depending  on  the 
exact  application.  Values  of  .02  to  .052  km  candles  were 
measured  in  fog  experiments  (Douglas  &  Booker,  p.  4- 
31).  For  pilot  landing  scenarios,  tests  yielded  values 
ranging  from  about  2  to  lOon  mile  candles. 

A  .scric.s  of  tests  was  run  at  Nantucket  Island  in  1941, 
comparing  transmissomctcroutput  with  visual  sighting.s. 
In  lhc.se  tests,  the  best  lit  was  obtained  not  by  using  a 
constant  illuminance  threshold,  but  by  using  an  illumi¬ 
nance  threshold  which  itself  is  a  function  ofvisibiliiy.  If 
S  is  the  threshold  obtained  at  1  unit  distance,  this  yields 
a  new  equation, 


s/v  =  ir/v^ 

(3.5) 

S  =  17'' /V 

(3.6) 

where  S  =  .052  for  visibility  in  km,  or  .084  for  visibility 
in  miles  when  I  =  25  candle  (Douglas  &  Booker,  p.  7- 1 7). 
This  relation  in  effect  states  tliat  although  the  illumi¬ 
nance  will  still  be  as  shown  in  Eq.  3.4,  the  effective 
illuminance  threshold  can  be  much  higher  in  low  visibil¬ 
ity  conditions  than  in  clear  air. 


Douglas  and  booker  iniiicaic  ihal  this  is  the  rclalion 
that  has  been  used  in  the  I  niicd  States  for  convertmg 
Iransniissometer  values  lo  night  visibility  since  lh;.i  time. 
'Fhe  tables  in  the  current  l  ederal  Meteorological  Hand¬ 
book  are  consistent  with  Equation  3.6. 

Thus  whereas  E(.|.  3.4,  Allard’s  law',  is  the  defining 
equation,  the  threshold  w  as  haind  in  at  least  one  scries  of 
tests  to  be  a  function  of  visibility.  'Hus  y  ielded  Eq.  3  6, 
which  is  the  equation  that  is  used  in  eonvening  iransmis- 
someter  values  to  visibilits  TItis  eijuation  w  ill  be  called 
the  transmissonieter  visibility  in  this  section. 

3.1.2  Comparison  with  Koschmictlcr  Visibility 

Equation  3.6  yields  a  somewbat  dilTcent  rel.ition 
between  visibility  and  transmittance  than  is  mlicivni  in 
the  daytime  definition  of  \  isibilily.  In  the  daytime,  the 
visibility  is  defined  as  the  distance  at  which  a  black  object 
of  certain  size  may  be  .seen  against  the  bori/on  sky .  1 1  can 
be  show-n  that  this  reduces  to  the  ctituilion 

r  7'  {3.?t 

when  thc.se  ideal  conditions  arc  mainlained.  In  Eq.  3.7, 
e  is  the  hum;m  contrast  threshold,  T  is  transmissivity,  and 
V  is  visibility.  This  is  another  fomi  of  Ko.schmieder's 
law  (Douglas  &  Booker  Seen  4.1.2).  The  WMO  recom¬ 
mends  a  threshold  of  .05.  Thi.s  is  the  rclalion  used 
currently  in  scatter  meters.  Based  on  historic  data, 
transmis.someters  u.se  Eq.  3.7  with  a  e  threshold  value  of 
.055. 

The  day  (Koschmieder)  and  night  (transmissometer) 
visibilities,  as  defined  by  Eq.  3.7  and  3.6,  arc  compared 
in  Figures  8  and  9.  Fugure  8  shows  the  visibility  values 
as  a  function  of  the  extinction  cocfncicni  (a  in  Eq.  3.3). 
F-igurc  9  show's  the  night  visibility  as  a  function  of  day 
visibility.  In  Eig.  8.  for  a  given  value  of  extinction 
cocITicicnl,  the  day  and  night  visibility  values  may  be 
quite  diflcrcni.  'I  bis  implies  that  the  value  of  visibility 
makes  an  abrupt  change  at  du.sk,  even  if  there  is  no 
change  m  the  atmosphere.  This  is  an  artifact  of  the 
current  approach  to  tire  definition  of  visibility;  rather 
than  being  a  tcmi  wiiich  is  a  well  defined  measure  of  the 
atmosphere,  it  is  a  tenn  which  attempts  to  predict  human 
response  .41  dusk,  as  the  available  targets  change,  and 
thchumanrcspon.se  changes,  the  visibility  determination 
changes. 

3.1.3  Common  Usage  of  the  Vi.sibilily  Equations 

Most  users  appearconsistcnl  in  using  the  Koschmieder 
equation  during  the  day.  During  the  day,  transmissom- 
eters define  visibility  in  tenns  of  tran.smissiviiyusing  Eq. 
3.7,  Uic  Koschmieder  cqutilion,  and  a  threshold  of  ,055 
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Fig.  8  Koschmieder  Visibility  (Eqn.  3.7)  and  Transmissometer  Night 
Visibility  (Eqn.  3.6)  as  a  Function  of  Extinction  Coefficient 
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Transmissometer  Night  Visibility  as  a  Function  of  Koschmiecter  Visibiiity 


Scatter  meters  use  a  threshold  of  ns  The  HSl  daytime 
algorithms  arc  also  consistent  with  lliis  dcHniiion  and  a 
threshold  of  .05. 

At  night,  even  though  Allard’s  law  is  a  defining 
equation,  both  Koschmteder  visibility  equation  (Eq.  3.7) 
and  transmissometer  visibility  equation  (Eq.  3.6)  arc 
used,  as  shown  in  Fig.  10.  Atnight.transmissomctcrsusc 
Eq.  3.6  with  a  threshold  of  .084  (mi).  The  Europeans  use 
Et].  3.7  and  a  threshold  of  .05  at  night.  Even  though  the 
transmis.soniclcr  output  is  apparently  irsed  only  in  low 
visibility  conditions,  tJic  original  test.s  indicate  Eq.  3.6 
should  also  be  a  rca.sonable  relation  to  use  in  clearer  air. 


Day  Visibility 

> 

l- 

II 

u 

Koschmieder 

£  =  .05 

WMO,  Scatter  Meters 

E=  .055 

Transmissometers 

Night  Visibility 

S  =  ITVV 

^  Allard 

S  =  ITVV 

Transmissometers 

t  =  T" 

European;  scatter  meters 

Fig.  10  Visibility  Equation  Usttgc 


The  weather  reports  use  mostly  visually  determined 
valuc.s,  which  arc  related  to  the  state  of  the  atmosphere, 
but  arc  certainly  not  a  tight  measure  of  the  state  of  the 
atmosphere.  Undcridcal conditions, the  psychophysically 
defined  visibility  may  be  theoretically  related  to  the  state 
of  the  atmosphere,  fn  practice,  however,  the  human  often 
uses  targets  that  arc  not  ideal.  If  a  non-ideal  target  at  a 
given  range  i  s  just  visi  bl  e,  then  that  range  i  s  defined  as  the 
visual  range  for  that  target.  That  visual  range  as  dis¬ 
cussed  by  Douglas  &  Booker,  p.  2-17,  is  not  the  same  as 
the  visibility. 

The  level  of  the  background  luminance  has  a  strong 
impact  on  the  ability  of  a  human  to  .see  a  given  light. 
However  tile  definition  of  night  visibility  may  be  consid¬ 
ered  to  include  the  stipulation  that  the  light  he  .seen 
against  an  c.sscnlially  black  background. 

The  requirement  for  black  background  is  not  staled 
explicitly,  however  Eq.  3.4  and  3.6  arc  only  valid  if  the 
background  is  c.sscnlially  black.  Since  this  is  the  normal 
dclming  equation,  and  is  also  u.scd  in  applying  iran.smis- 
someter  data,  it  may  be  considered  part  of  the  normal 
definition. 


Douglas  and  Bookcrtfi  4  -  17|discussihc  implications 
of  a  non-black  background  ui  greater  detail.  In  some 
visual  problems,  such  as  runway  approach  hgiits,  it'.c 
background  is  less  than  Iff  of  the  source,  in  others,  such 
as  with  taxiway  lights,  the  txickground  can  be  significant 
It  would  be  feasible  w  ith  the  I  LSI  to  correct  Eq.  3.6  ionlie 
background  that  obtains  at  the  lime  o!  rneasurenteni 
(Douglas  &  Hooker  p  4-17),  n  order  to  predict  l)Ow  a 
human  might  see  25  candle  lights  against  inis  parncuiar 
btickground. 

We  propose  that  llic  llSl  compute  the  Ko.schmieder 
visibility  tlai.  3.7  with  a  thre.shold  of  .05)  as  the  best 
measure  of  the  state  of  the  atmosphere,  in  addition,  the 
night  transmissometer  visibility  definition  {Eq.  .^.6  with 
a  threshold  of  .084  {mip  could  be  used  to  provide  an 
optional  output. 

Considering  the  multiplicity  of  defining  expressions 
and  icnns,  as  noted  in  the  preceding  paragraphs,  it  should 
be  noted  that  the  terms  Runway  Visibility  Value  (RVV), 
Runway  Visual  Range  (RVR),and  Meteorological  Opti¬ 
cal  Range  (MOR)  arc  other  parameters  closely  related  to 
visibility.  A  full  discu.ssion  of  these  concepts  i.s  beyond 
the  scope  of  thi.s  report.  A  variety  of  definitions,  includ¬ 
ing  psychophy.sical  definitions,  arc  given  in  Douglas  and 
Booker  (1977)  and  the  Federal  Meteorological  lland- 
iHiokNo,  1  (1988). 

3.2  Measuring  the  Optical  Attenuation  at  Night 

The  above  defining  equations  relate  the  visibility  to 
ihciransmissiviiyofihc  atmosphere.  This  transmissivity 
may  be  determined  with  the  HSl  u.sing  cilherof  two  main 
approaches:  basing  the  computations  on  absolute  mea¬ 
surements,  and  basing  them  on  relative  measurements 
Gordon  (1989)  suggests  using  measurements  of  the  rela¬ 
tive  illuminance  (apparent)  of  lights  of  equal  inicmsily 
(which  is  an  inherent  parameter).  Theoretically,  it  should 
be  feasible  to  determine  the  extinction  from  any  .set  of 
lights  of  known  inherent  relative  output.  This  section 
di.scus.sc.s  first  the  absolute  approach,  ilicn  the  relative 
approach. 

The  equations  will  differ  somewhat  depending  on 
whether  the  targets  are  sub-pixel  in  si/c,  nr  larger.  If  the 
targets  arc  largcrthan  a  pixel,  then  they  may  be  treated  as 
an  extended  source,  and  one  uses  the  measured  ab.solutc 
or  relative  luminance.  If  the  targets  arc  sub-pixel,  they 
act  as  a  point  source.  In  this  ca.se  the  luminance  is 
undefined  and  one  uses  the  intensity.  Section  3  develops 
the  equations  for  llio  extended  source. 

Note  that  in  lliis  section  either  the  term  radiance  or 
luminance  could  be  used.  These  terms  arc  noicquivaleni. 
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however  Oic  beam  transmittance  may  lx  determined 
from  either  parameter. 

3.2,1  Ab.solutc  Measurement  Approach 

With  an  extended  source,  the  absolute  approach  to 
determining  transmissivity  is  theorelically  quite  simple. 
The  apparent  luminance  (Ly  of  a  light  at  known  distance 
is  measured.  If  the  inherent  luminance  (Lo  is  known  (or 
derived  from  clear  night  measurements),  then  the  beam 
transmittance  may  be  found  from 


3.2.2  Relative  .Measureinent  .Approach 

In  this  sect  ion,  wederi'.c  the  ecjuaiiotis  which  could  lx 
used  todenve  visibility  from  relative  measurements.  As 
in  Section  3,2.1.  we  will  assume  extended  sources;  the 
changc.s’  to  point  .sources  are  given  in  Section  3. 2. .3, 

Consider  a  set  ol  sources  with  known  relative  inherem 
illuminance.  We  will  use  "L''  for  absolute  illuminance 
and“N”  for  relative  illuminance.  Then  source  “i”  has  an 
absolute  inherent  illuminance  given  by 


T,  =  ,L,/,L„  .  (3.8) 

The  transmissivity  may  be  derived  from  tlic  measured 
luminance  by 

r  =  T/^  =  (,U/.L„)'^' ■  (3.9) 

The  Koschmieder  visibility  equation  3.7  may  be  ex¬ 
pressed  as 

V  =  Ine/InT  .  (3.10) 

Substituting  Eq.  3.9  for  T,  the  equation  for  determining 
visibility  given  the  ab.solutc  luminance  values  for  a  target 
becomes: 

y  =  lne/ln(,L,/,L„)  .  (3.11) 

To  determine  night  tran.smi-s.somctcr  vi.sibilily  from 
the  ab.solutc  measurements,  one  may  use  Eq.  3.6,  which 
gives  the  night  transmis.somctcr  visibility  as  a  function  of 
ihctraasmi.ssiviiy.  The  iransmis.sivity  is  derived  from  the 
measured  values  using  Eq.  3.9,  Once  this  traasmissivity 
is  determined,  the  night  transmissomclcr  visibility  is 
determined  iteratively  from  Eq.  3.6. 

It  should  be  noted  that  the  night  visibility  is  defined 
only  in  terms  of  point  sources  (Eq.  3.6  implies  a  point 
source).  It  might  .seem  odd  to  measure  it  using  extended 
sources.  The  rca.son  this  is  possible  is  that  the  transmit¬ 
tance  may  be  determined  using  cither  point  sources  or 
extended  .sources.  Equation  3.6  then  gives  an  approxi¬ 
mate  indication  of  how  far  a  human  could  see  a  point 
source,  given  an  atmosphere  with  the  given  transmissiv¬ 
ity,  (The  distance  over  which  a  human  could  see  the 
extended  source  is  not  Uic  same  as  the  night  transmi.ssom- 
clcr  vi.sibility  defined  by  Eq.  3.6.) 

Whereas  the  ab.solutc  approach  is  theoretically  quite 
convenient,  it  requires  a  very  stable  system,  and/or  a 
mcansof  correcting  forsystem  absolute  response  changes. 
For  this  rca.son,  it  may  be  preferable  to  use  a  visibility 
determination  which  docs  not  depend  on  stability  of  the 
sensor  over  long  periods. 


.L,  =  C;  ,\„  (.3,12) 

where  C’l  is  an  unknown  cnnsiam  that  is  the  same  lor  all 
of  tlie  sources.  We  are  able  to  measure  ihe  relative 
apparent  illuminance  given  by 

N  =  C,  L  (3.13) 

where  C2  is  an  unknown  constant  that  may  vary'  from 
image  to  image,  as  the  camera  absolute  response  changes. 
The  apparent  illuminance  is  related  to  the  inherent  illu¬ 
minance  by 

,L,  =  7'',L„  (3.14) 

and  therefore  the  apparent  relative  illuminance  is  related 
to  the  inherent  relative  illuminance  by 

.N,  =C2.L,  =Cj(r  ,l.„)  =  C,C:r  .  (3.1.S) 


If  one  mca.sures  tlic  apparent  relative  illuminance  of 
two  sources,  with  known  inherent  relative  illuminance, 
one  may  ratio  them  to  derive  the  iransmi.s.sivity,  as 
follows.  First,  the  ratio  ol  the  measured  relative  values  is 
given  by 

,N,  _  C.C^  r'  ,N,  ,N„ 

2N,  C,C2/'’2N„  jNo  (3  16) 


Where  both  targets  arc  from  the  .same  image  so  that  C2  is 
constant.  Taking  the  logs  of  both  sides  yields 


In 


2N, 


=  (r,  -r2)  In/'  +  ln 


iis 


(.3.17) 
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(3.18) 


Eq.  3. 1 8  may  be  substituted  directly  into  the  equation 
for  Koschmieder  visibility,  as  exprc.sscd  in  Eq.  3.10. 
This  yields  an  equation  for  deriving  visibility  from  rela¬ 
tive  luminance  values: 
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V  =  (r, -fj)  Ine/ln 


|Nr  2^0 
aN,  ,N„ 


(3.19) 


For  night  transmissometcr  visibility,  wc  derive  the 
transmissivity  from  Eq.  3.18 

i/(r,-,2) 


T  = 


2N,  ,N„J 


(3.20) 


and  then  use  Eq.  3.6  to  iteratively  detennine  the  night 
visibility.  As  discussed  in  the  previous  section,  this  is  the 
night  visibility  of  point  sources,  even  though  it  has  been 
derived  indirectly  by  using  transmissivity  determined 
from  measurements  of  extended  sources. 

3.2.3  Summary  of  the  Visibility  Equations 

As  shown  in  Fig.  1 1 ,  there  arc  several  choices  for 
defining  equations  to  use  with  the  night  WSI.  First,  there 
is  the  Koschmiedcr’s  law,  used  during  daylight  and  also 
used  in  many  cases  at  night.  Next  is  Allard’s  law,  from 
which  the  night  transmissometer  visibility  definition 
arises.  We  recommend  using  the  Koschmieder  visibility 
with  the  HSI,  with  the  night  transmissomctervisibility  as 
an  option  to  be  added  at  a  later  time. 

On  the  right  side  of  Fig.  1 1  are  shown  the  methods  of 
using  the  HSI  to  determine  the  transmissivity  which  is 
input  to  those  defining  equations  on  the  left  half  of  the 
figure.  The  top  equation  is  that  used  with  the  HSI  during 


the  day,  and  is  based  on  the  measured  contrast.  At  night, 
the  transmissivity  may  be  determined  from  either  the 
absolute  or  relative  radiance  of  lights,  as  di.scussed  in 
.sections  3.2.1  and  3,2.2,  yielding  the  remaining  equa¬ 
tions  shown  in  the  figure.  Wc  feel  that  the  relative 
approach  should  be  the  most  accurate  detciTnination  (due 
to  camera  long-term  drift),  but  feel  both  approaches 
should  be  tested.  Thus  the  recommended  equations  are 
Eq.  3.19  for  the  primary  ctjuation,  and  Eq.  3.1 1  as  the 
alternative  lest  equation. 

3.2.4  Conversion  to  Point  Sources 

If  a  given  source  is  smaller  than  a  pixel,  then  it  is 
effectively  a  point  source.  In  this  ease,  the  luminance  of 
the  source  is  undefined,  and  the  loss  of  flux  is  measured 
in  terms  of  illuminance.  The  beam  transmittance  is 
unchanged,  but  whereas  the  loss  of  luminance  is  given  by 
Eq.  3.8,  the  loss  of  illuminance  is  given  by  Eq.  3.2. 

One  can  show  that  the  equations  given  in  Section  3.2 
may  be  convened  from  extended  to  point  .sources  by 
using  the  apparent  illuminance,  tEr,  in  place  of  the 
apparent  luminance  (Nr,  and  by  using  (1  /  in  pi  acc  of  the 
inherent  luminance  iNo.  In  the  above,  the  intensity  is 
denoted  il  rather  than  iIq.  since  intensity  is  by  definition 
an  inherent  property.  Thus,  Eq.  3.19  becomes 


V  =  (r,  -r^)  Inc /In 


(3.21) 


Here  E  and  I  should  be  understood  to  be  relative  values 
of  illuminance  and  intensity  and  Eq.  3.1 1  becomes 


V  =  fine /In 


(3.22) 


In  practice,  if  the  range  to  a  target  will  not  be  changing, 
it  is  probably  most  convenient  to  use  equations  analo¬ 
gous  to  Eq.  3. 19  and  3. 1 1  where  the  apparent  property  is 
the  apparent  illuminance,  and  the  inherent  property  i.s  tftc 
illuminance  that  would  be  observed  at  the  target  range 
given  zero  attenuation.  In  this  way,  the  system  can  be 
programmed  with  tire  .same  equations,  whciltcr  die  sources 
arc  point  or  extended,  as  long  as  care  is  taken  in  under¬ 
standing  the  appropriate  inherent  property  to  input. 

3.3  Hardware  Adaptations  for  Visibility  at  Night 

In  order  to  dclcmiinc  visibility  at  night,  it  was  neces¬ 
sary  to  have  a  system  with  sufficient  sensitivity  for 
imaging  lighLs  of  opportunity  at  night.  As  di.^cu.s.scd  in 
Johnson  cr 0/(1990),  initial  tests  with  an  intensified  CID 
system  did  not  yield  adequate  stability  and  image  quality. 
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It  was  decided  to  utilize  the  injeciioii  inhibit  capability  of 
the  CID  cameras,  in  which  tlie  signal  is  allowed  to 
accumulate  formultiple  frameperiods  before  rcadingout 
the  accumulated  signal  (before  injecting  the  electrons 
into  the  substrate  for  readout).  This  necessitated  three 
changes:  changing  the  camera  usage,  changing  to  a  new 
imaging  board,  and  adding  an  interface  circuit. 

3.3.1  CID  Injection  Inhibit 

The  following  description  of  the  CID  injection  inhibit 
is  extracted  from  Johnson,  et  al,  (1990). 

The  sensorsub-system  u.scd  in  the  HSI  is  the  CID  27 1 0 
which  incorporates  CIDTEC’s  exclusive  Charge  Injec¬ 
tion  Device  (CID)  sensor  into  a  high-quality,  solid  state, 
monochrome  RS-17()  camera.  (Ref.  Vol.  II) 

“The  Charge  Injection  Device  is  an  X-Y  addressable 
image  sensor  organized  so  that  photon-generated  charge 
can  be  read  at  each  sensing  site  (pixel)  through  local 
charge  transfer.  Injection  into  the  underlying  substrate  is 
used  to  clear  each  pixel  of  signal  charge  and  start  a  new 
signal  integration  period  (television  frame  time).  These 
images  arc  fabricated  in  an  epitaxial  layer;  the  epitaxial 
junction  is  reverse  biased  and  used  as  acommoncollcctor 
to  remove  the  injected  charge  from  the  device.  The 
coascquences  uf  this  organization  arc  numerous.  Since 
charge  is  scased  locally,  there  arc  no  paths  along  which 
optical  overloads  or  defects  (black  or  white  faults)  can 
propagate.  Charge  transfer  losses  arc  not  cumulative. 
Pixels  can  be  contiguous,  leading  to  high  modulation 
transfer  function  and  high  quantum  efficiency.  The 
sensor  has  inherently  low  dark  current;  the  depletion 
region  of  each  pixel  is  smaller  than  the  photo  sensitive 
region  and  inversion  charge  can  be  used  to  quench 
surface  leakage  current.  Read-out  can  be  non-destruc¬ 
tive  allowing  on-chip  signal  processing,  and  the  X-Y 
addressable  array  can  be  random  accessed,”  (Grafingcr 
&  Michon) 

The  HSI  system  takes  advantage  of  the  inherent  char¬ 
acteristics  of  the  CID  2710  camera  system  to  address 
several  low  light  level  measurement  scenarios  without 
the  need  for  external  image  intensification.  This  is 
accomplished  by  implementing  the  Inject  Inhibit  func¬ 
tion  available  in  the  standard  camera  configuration.  In 
this  mode,  normal  destructive  readout  is  interrupted  by 
the  Inject  Inhibit  function  which  permits  the  imager  to 
integrate  for  longer  than  one  frame  time.  The  2710 
continues  to  integrate  as  long  as  the  Inject  Inhibit  func¬ 
tion  is  engaged,  facilitating  user  controlled  time  expo¬ 
sure.  Several  mock-up  circuit  boards  were  built  by  MPL 
to  enable  the  selection  of  integration  intervals  ranging 
between  2  and  2(X)  frames.  At  standard  video  rates  of  30 


frames  per  ncco-kI,  ituegraiion  intcivals  represent 
roughly  60  milliseconds  upward  to  6  seconds, 

3.3.2  The  VS-  !(XI  Image  Board 

Under  integration,  Uic  CID  camera  outputs  a  signal 
which  remains  in  RS170  fonnat,  however  during  tlie 
integration  interval,  the  output  frames  on  the  K.S17() 
signal  have  0  analog  signal,  i.e.  are  black.  I-igure  12 
shows  a  simplified  rendition  of  the  R.S170  signal  at 
vertical  frame  rate.  Al  0  integration,  i.e.  normal  video 
rate,  the  ClDouiputs  each  field  separated  by  venical  sync 
information,  with  two  fields  creating  a  full  frame.  (Each 
field  consists  of  the  output  from  either  even  or  odd  scan 
lines,  with  horizontal  sync  infomialion  between  lines,  ref 
Varah  and  Sprink.  1991 .)  With  integration  set  to  1.  the 
camera  integrates  during  an  extra  frame,  during  which 
time  the  active  video  output  is  0.  This  is  followed  by  a 
frame  with  the  integrated  video  signal  read  out.  Simi¬ 
larly,  as  shown  in  Fig.  12,  at  higher  integration  settings, 
the  video  output  is  0  during  integration,  followed  by  a 
read-out  frame. 

The  signal  is  read  by  tlie  control  system  through  an 
imaging  board,  which  pcrfomis  the  following  on-board 
functions:  A/D  conversion,  frame  memory,  image  pro¬ 
cessing,  and  D/A  fordisplay  on  the  monitor.  The  daytime 
HSI  and  WSI  both  u.sc  an  Imaging  Technology  FG-100 
board.  For  night  operation  with  the  injection  inhibit 
function,  it  was  necessary  to  change  to  the  VS- 1 00  board . 
The  primary  new  feature  of  this  board  u.sed  in  this 
application  is  the  external  trigger.  That  is,  this  board 
allows  use  of  an  external  trigger  for  frame  grab,  so  that 
the  readout  frame  could  be  grabbed,  rather  than  one  of  the 
black  frames  output  during  integration. 

3.3.3  The  Integration  Control  Circuit 

In  order  to  fully  synchronize  the  VS-100  board  with 
the  camera  injection  inhibit,  it  was  necessary  to  design 
and  build  an  integration  control  circuit.  This  circuit  is 
documented  in  Varah  and  Sprink  ( 1991). 

The  control  interface  allows  control  of  the  integration 
period  cither  through  computer,  or  manually  via  thumb 
wheels.  The  control  circuit,  synced  by  the  vcnical  sync 
from  the  camera,  .sends  a  signal  to  the  camera  to  start 
integration.  At  the  same  time,  it  begins  counting  integra¬ 
tion  periods.  On  completion  of  the  desired  integration 
interval,  it  sends  a  signal  to  the  frame  board  to  grab  the 
next  output  frame,  and  sends  a  signal  to  the  camera  to  stop 
integration.  The  above  is  a  conceptual  overview;  in 
practice,  the  control  is  considerably  complicated  by  the 
need  for  appropriate  delays  for  coordinating  the  camera 
and  frame  board. 


11 


odd  even 

field  field 


full  frame 


vertical  sync 

I  J 

\r 


IT 


LT 


integration  readout 


integration  readout 


I  *'  ■  I  I  i-i  III  p 


integration 


readout 


If 


integration 


readout 


Inj  Inh 
Setting 


Fig.  12  Conceptual  illustration  of  output  RS170  Video  Signal 
(Vertical  frame  rate)  at  various  injection  inhibit  sellings 


3.3.4  The  Auto-iris  Inhibit 

Normally  the  auto-iris  reads  the  output  video  signal, 
and  adjusts  the  aperture  opening  to  keep  the  signal  at 
some  preset  level,  nominally  mid-scale.  There  were  two 
concerns  with  the  use  of  the  auto  iris  at  night.  First,  if  the 
camera  integrated  long  enough  for  the  signal  to  become 
reasonably  bright,  the  iris  might  begin  to  close  in  re¬ 
sponse.  This  would  be  counter-productive;  it  makes  no 
sense  lo  use  integration  for  incrca.scd  signal,  then  de¬ 
crease  the  signal  with  an  aperture.  Secondly,  depending 
on  the  integrati  on  period  used  with  the  i  ri  s ,  the  re  m  i  ght  be 
unpredictable  response  tc  the  changing  video  output 
during  integration,  as  the  output  changes  from  black 
during  integration  to  grey  during  readout  (as  .shown  in 
Fig  12). 

In  order  to  defeat  lhc.se  problems,  it  was  decided  to 
control  the  auto  iris  such  that  any  lime  injection  inhibit  is 
used,  the  auto  iris  is  disabled  and  forced  fully  open. 
When  an  injection  inhibit  setting  of  0,  or  video  rate  is 
used,  the  auto  iris  is  enabled. 

All  of  the  above  hardware  changes  have  been  built  and 
tested.  They  perform  normally,  and  arc  considcicd 
operational  at  this  lime. 


3.4  Software  Adaptations  for  Visibility  at  Night 

There  were  very  significant  changes  to  the  HSIAVSl 
software  during  the  last  year.  As  discussed  in  section 
3.4.1,  software  to  enable  night  data  acquisition  was 
written.  Section  3.4.2  discusses  the  changes  to  the 
daytime  code,  for  compatibility  with  the  new  VS-1(X) 
imaging  board.  Additional  changes  to  die  daytime  code, 
to  enhance  its  capabilities,  arc  discussed  in  Section  4, 
whichdi.scu.sscs  refinements  to  the  HSI  daytimccapabi!  i- 
lics. 

3.4, 1  Night  HSI  Programs 

Three  programs  were  written  for  handling  data  acqui¬ 
sition  with  integration.  The  primary'  program  for  night 
data  acquisition  is  Program  GETNITF..  This  program  is 
dc.signcd  lo  run  unattended,  lo  let  the  user  grab  night 
imagery  and  store  the  images  on  Exabyte  tape.  The  user 
maycnicrslarl  time  and  stop  time,  and  acquire  images  at 
a  variety  of  integration  levels  automatically. 

Program  INTGR8  is  a  utility  designed  for  u.sc  in 
calibration  and  oilier  camera  tests.  This  program  allows 
the  uscrlo  enter  the  integration  level  cither  via  keyboard 
ormanually,  and  acciuire  and  .save  imagery.  The  integra¬ 
tion  level  may  be  varied  during  the  program, 
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Finally.  Program  VS  UTIL  is  a  general  udlily  program 
for  accessing  ihc  VS  100  board,  1:  was  quite  preliminary, 
however  it  may  still  be  used  for  accessing  the  data  stored 
on  tape,  when  a  machine  with  a  V^SIOO  hoard  in  place  is 
u.sed.  This  is  not  normally  used  in-house,  as  imagery 
saved  with  any  of  the  above  programs  may  tx;  acces.scd 
on  any  of  the  in-house  computers  with  FG-lOO  boanls. 

3.4.2  Adaptation  of  the  Day  HSI/WSl  programs  for 
VS-lOO 

Changing  from  the  FG 100  image  board  to  the  VS-  1(X) 
board,  to  enable  using  injection  inhibit  for  on-betard 
integration  with  the  CID,  had  one  unfortunate  result:  tlie 
software  had  to  be  changed,  even  for  daytime  use,  be¬ 
cause  the  VS- 100  board  is  not  downwardly  compatible 
with  the  FG-100  board.  As  a  result,  all  of  the  daytime 
code,  for  both  WSI  and  HSl,  had  to  be  changed.  These 
programs  are  listed  below. 

Primary  Programs 

Setupvis  -  Creates  azimuth  selection  flic,  allows  user 
to  change  input  parameters  (range,  thresh¬ 
old  and  inherent  contrast), 

VistcxlS  -  Automated  visibility  determination 
program. 

Wsiten  -  Cloud  cover  determination  program. 


3.5  CII)  Camera  llehavior  under  Inlegrati«/i 

In  order  to  acquire  tlic  night  data,  it  was  necessary  to 
use  the  CID  camera  in  integration  mode,  as  discussed  in 
Section  3.3.  A  sericsof  tests  and  radiometric  calibrations 
were  run  irn  the  CID  cameras  to  deienninc  whether  this 
mode  of  operation  introduced  any  signal  peculiarities 
which  must  be  accounted  for  in  the  data  handling. 

Several  in-house  camer.is  were  tested,  and  found  to  be 
deficient  in  a  variety  of  ways,  even  though  they  behaved 
well  under  normal  video  mode  use.  Since  the  faults 
appeared  to  bo  related  to  aging,  two  new  cameras  wore 
purchased.  These  were  found  to  behave  adet|uatcly.  This 
section  first  describes  the  I'lchavior  of  the  new  sensors, 
and  then  discusses  the  types  of  faults  which  occurred  in 
those  cameras  which  had  begun  to  age. 

The  CID  camera  is  reasonably  linear,  in  the  sense  that 
a  fractional  change  in  input  flux  yields  a  similar  frac¬ 
tional  change  in  output  signal.  Fig.  13  shows  the  results 
of  a  linearity  calibration,  acquired  at  normal  video  rate. 
The  system  is  quite  linear  over  most  of  the  span,  with 
some  decrease  in  signal  relative  to  the  expected  signal  at 
the  bright  end.  This  is  quite  typical  of  the  CID  .scasors. 
If  the  sensor  is  used  in  the  visibility  system  without 
correcting  for  this  slight  non-linearity,  the  rc.sulting  error 
is  typically  7%,  when  the  horizon  radiant  signal  is  about 
200.  As  mentioned  earlier,  wc  recommend  changing  the 
current  HSl  algorithms  by  applying  the  linearity  calibra¬ 
tion  results  to  the  imagery  to  correct  for  the  non-linearity. 


Fixhomc 

-  Fine  tunes  home  position. 

AV5 

-  Program  is  now  called  Snglvis.  Computes 
visibility  for  a  range  of  inherent  contrasts. 

Scevis 

-  Displays  image  di.sk  file  to  RGB  screen. 

Savimg 

-  Allows  user  to  digitize  an  image,  then 
save  it  to  disk. 

Savrgb 

-  Allows  user  to  save  image  currently  on 
RGB  screen  to  di.sk. 

Vicwvis 

-  8mm  data  tape  created  by  Vistex  15. 

Programs  such  as  Move  which  move  the  rotary  tabic 
to  a  user  .selected  azimuth  and  Printvis  which  print-s  and 
translates  visibility  data  files  into  a  spreadsheet  fomial 
arc  not  affected  by  the  VS  100  board. 

The  changes  to  the  daytime  programs  to  enhance  to 
enhance  their  capability  are  discussed  in  Section  4. 


When  the  sensor  is  calibrated  using  tnicgralion,  an 
abnormal  feature  becomes  immediately  apparent  onmany 
of  the  cameras,  including  one  of  the  two  new  ones.  TTie 
alternate  field  lines  have  slightly  different  magnitude; 
that  is,  the  odd  field  and  even  field  have  slightly  different 
signals  (ref  Fig.  12).  Whereas  it  is  part  of  our  normal 
procedure  to  adjust  the  alternate  field  lines  at  normal 
video  rate  in  the  electronic  calibration,  there  is  not  a  way 
to  adjust  this  under  integration.  A.s  an  example,  the  new 
camera  with  Uiis  problem  had  a  standard  deviation  in  the 
signal  over  a  20  x  20  pixel  square  of  2  counts  in  video 
mode,  and  about  10  counts  at  integration  level  I. 

As  shown  in  Fig.  12,  llic  output  signal  is  nominally 
zero  during  the  integration  period.  With  tho.se  sensors 
.showing  the  alternate  field  offset,  there  is  a  non-zern 
signal  for  one  of  llic  fields  even  during  the  inicgration 
period.  For  this  rca.son,  wc  would  recommend  further 
cvalualionof  potential  wavs  to  corrccUhis  ahcmaic  field 
problem  prior  to  purchase  of  any  additional  units.  It  may 
be  preferable  at  this  point  to  convert  to  a  CCD  .sensor,  if 
one  with  better  operating  characteristics  is  available  at  a 
similar  price. 
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Fig.  1 3  Results  of  the  Linearity  Caiibration,  showing  the  silght 
non-linearity  of  the  CIO  Sensor  at  0  Integatlon 


Results  from  asamplc  calibration  of  signal  vs.  integra¬ 
tion  period  are  illustrated  in  Fig.  14.  The  pixels  continue 
to  receive  flux  during  the  readout  interval,  as  well  as 
during  the  integration  interval,  so  the  expected  signal 
should  be  proportional  to  (1  + 1 ),  where  I  is  the  integration 
interval.  In  Fig.  14,  the  straight  line  represents  this 
expected  relationship.  The  measured  rcspon.se  is  quite 
close  to  this  expected  response.  The  lower  curve  in  Fig. 
14  will  be  discu.sscd  later  in  this  .section. 

Measurements  similar  to  those  shown  in  Fig.  14  were 
also  acquired  for  a  variety  of  other  signal  levels  and 
integration  periods.  It  was  determined  that  whereas  this 
camera  is  well-behaved  out  to  an  integration  of  10 
periods,  the  signal  is  non-linear  beyond  an  integration 
level  of  20  periods,  even  when  the  signal  is  low  enough 
to  be  well  onscalc.  Beyond  20,  the  signal  becomes  quite 
non-linear,  with  an  output  lower  than  expected  based  on 
the  (I  +  1)  relationship. 

Thus  the  new  CID  sensors  were  reasonably  linear  in 
terms  of  signal  vs.  input  flux,  and  quite  linear  in  terms  of 
signal  vs.  integration  period  out  to  10  periods  integration. 
The  CID  performance  in  other  respects  was  quite  ad¬ 
equate,  except  that  .some  new  .sensors  apparently  have  a 
field  offset  in  integration  mode  which  degrades  the 
quality  of  the  measurements. 

The  response  of  the  CID  scn.sors  to  aging  is  also  a 
serious  concern  with  respect  to  their  use  under  integra¬ 
tion.  Wc  had  in-house  several  sensors  of  varying  ages, 
which  were  still  quite  good  for  normal  full  video  rate 
usage.  That  is,  the  dark  levels  were  reasonably  low, 
standard  deviations  low,  and  other  features  such  as  ramp¬ 
ing  and  field  offset  low.  However,  under  integration, 
none  of  these  older  sensors  (same  model,  longer  usage) 
performed  adequately.  One  camera  had  full  dark  levels 
which  became  quite  high  under  integration.  Another  had 
a  severe  field  offset  problem. 

The  most  common  problem  with  these  older  CID 
sensors  was  that  the  behavior  with  integration  was  no 
longer  linear.  The  lower  curve  in  Fig.  14  illustrates  the 
rolloff  feature  vs.  integration  period,  which  was  common 
to  the  older  cameras.  Not  only  was  this  rolloff  quite 
camera  dependent,  but  it  was  scene  dependent.  The 
example  given  was  acquired  with  a  uniform  image,  i.c. 
one  brightness  level  over  the  full  array.  The  rolloff  was 
much  more  extreme  when  only  one  portion  of  the  array 
was  illuminated.  Thus  the  calibration  becomes  both 
scene  and  integration  level  dependent  in  the  older  sen¬ 
sors. 

Wc  do  not  know  what  aging  rate  is  critical  for  the 
sensors.  Wc  feel  that  the  problem  may  be  due  to  recom¬ 


bination  of  tlic  clcciron-holc  pairs  at  the  sites  of  crystal 
defects,  and  that  a  measure  of  the  background  dark  level 
may  be  in  reasonable  indication  of  this  aging  process  in 
the  crystal.  However  this  is  .surmise  only  at  this  point 

It  is  fair  to  say  that  the  Cf  D  sensors  arc  quite  effective 
for  daytime  visibility,  and  arc  adequate  for  night  tests 
required  to  further  develop  the  concept  of  dcicmtining 
visibility  at  night.  Howeverthese  sensors  arc  not  (^piim.il 
over  the  long  term  for  night  visibility  acquisition,  due  to 
theirunfavorable  response  under  integration  after  aging. 
Slow  scan  CCD  sensors  have  more  than  adequate  sensi¬ 
tivity  for  both  night  and  day  visibility  applications,  at  a 
higher  cost.  A  high  priority  should  be  evaluation  of  iltc 
relative  cost  and  merits  of  the  variety  of  .solid  state 
.sensors  currently  available  for  this  task. 

3.6  Sample  Night  Visibility  Computations 

To  provide  a  sample  computation  of  the  techniques 
discussed  in  the  previous  sections,  wc  acquired  data 
during  the  daylight  and  evening  hours  on  a  clear  day  and 
a  foggy  day.  First,  the  sample  daytime  images  extracted 
for  study  were  used  to  derive  daytime  visibility  in  the 
normal  manner.  The  results  were  about  26  miles  on  the 
cicarday,  and  2  miles  on  the  foggy  day.  as  shown  in  Tabic 
1. 

Table  1 

Day  Visibility  Sample  Values 


Date 

Time 

Azimuth 

Visibility 

16  Jan 

234224 

80 

28  mi 

16  Jan 

235651 

110 

25  mi 

23  Jan 

001224 

80 

<  2.6  mi 

23  Jan 

002651 

110 

2.0  mi 

In  order  to  evaluate  the  night  data,  we  first  had  to 
estimate  the  inherent  signal,  i.c.  the  signal  corresponding 
to  1CX)%  transmittance.  To  do  this,  we  assumed  that  at 
night  on  the  16lh  the  visibility  was  still  about  26  miles 
and  used  this  data  to  determine  the  inherent  signal  from 
selected  lights.  Then  using  this  approximate  inherent 
value,  wc  computed  the  visibility  for  23  Jan  at  night, 
using  the  measured  night  signals. 

Wc  were  able  to  find  4  targets  which  were  onscalc 
both  nights  at  the  same  integration  period,  at  azimuth  80. 
At  azimuth  1 10.  there  were  no  targets  onscalc  both  nights 
with  the  same  integration  period,  but  the  use  of  different 
integration  periods  allowed  selection  of  two  targets  with 
onscalc  data  on  both  nights. 
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Results  of  the  Calibration  showing  camera  response 
to  Increasing  the  Injection  Inhibit  setting 
(Increasing  the  integration  period) 


The  compulation  of  the  inherent  signal  is  shown  in 
Table  2.  In  Table  2,  data  have  been  extracted  from  the 
clear  night  images.  The  background  signal,  a  combina¬ 
tion  of  path  radiance  and  .sensordark  level,  wascxtracted 
a  few  pixels  below  the  light,  and  liic  signal  corrected  for 
the  background  is  computed.  Given  the  range  to  the 
target,  the  approximate  transmittance  may  be  estimated, 
assuming  a  visibility  of  26  miles  (the  visibility  measured 
earlier  during  the  day).  From  this  computed  transmit¬ 
tance  and  the  background  corrected  signal,  we  may 
compute  the  signal  which  should  occur  if  the  beam 
transmittance  were  100%.  This  is  labeled  the  “pseudo 
inherent  signal”.  This  label  reflects  the  fact  that  since 
many  of  these  lights  may  be  point  sources,  the  p.scudo 
inherent  signal  docs  not  refer  to  the  signal  at  zero  range; 


rather  it  refers  in  this  case  to  the  signal  at  the  given  range 
which  should  occur  if  the  transmittance  were  KX)%. 

Once  a  pseudo  inherent  signal  is  computed  for  each 
target,  the  data  from  other  nights  may  be  used  to  estimate 
the  visibility.  'I'his  procedure  is  illustrated  in  Table  3 
using  the  foggy  night  imagc.s.  In  this  table,  the  first  signal 
shown  is  the  pseudo  inherent  signal  from  Table  2.  Next 
is  the  raw  signal  and  adjacent  background,  which  is 
subtracted  to  yield  the  corrected  signal.  At  the  110 
azimuth,  the  corrected  signal  also  reflects  an  adju.stment 
for  the  integration  period  of  7,  to  correct  it  to  the  signal 
appropriate  for  an  integration  period  of  4.  The  visibility 
is  then  computed  directly  from  Eq.  3. 11,  the  derived 
visibility  equation  using  the  absolute  approach. 


Table  2 

Computation  of  Inherent  Signal  From  Clear  Night  Data  (16  Jan) 


Date 

Tgt# 

Raw  Sig 

Bkgnd 

Corr  Sig 

Range 

(mi) 

Est 

Trans. 

P.seudo 
Inh.  Sig. 

Azimuth  80,  Int  level  4,  Time  04 

0601 

16  Jan 

n 

174 

36 

138 

1.2 

.87 

159 

16  Jan 

Wm 

107 

29 

78 

1.1 

.88 

89 

16  Jan 

95 

28 

67 

1.1 

.88 

76 

16  Jan 

n 

217 

30 

187 

1.1 

.88 

212 

Azimuth  1 10,  Int  level  4,  Time  035019 

16  Jan 

5 

112 

28 

84 

2.8 

.72 

117 

16  Jan 

6 

182 

20 

162 

2.8 

.72 

225 

Table  3 

Computation  of  Foggy  Night  Visibility  (23  Jan)  Using  Absolute  Equation 


Date 

Tgt# 

Pseudo 
Inh.  Sig. 

Raw 

Signal 

Bkgnd 

Corr 

Sig 

Vis  (mi) 

Azimuth  80,  Int  level  4,  Time  040601 

23  Jan 

1 

159 

71 

28 

43 

2.8 

23  Jan 

B 

89 

52 

25 

27 

2.8 

23  Jan 

B 

76 

49 

25 

24 

2.9 

23  Jan 

B 

212 

39 

26 

13 

1.2 

Azimuth  1 10,  Int  level  4,  Time  042218 

23  Jan 

5 

117 

53 

28 

15.6* 

4.2 

23  Jan 

6 

225 

50 

30 

12.5* 

2.9 

♦Corrected  for  both  background  signal  and  Integration  7  vs  Integration  4  calibration 
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In  Tabic  3,  most  of  the  values  arc  reasonably  consis¬ 
tent  indicating  a  visibility  of  about  3  miles.  This  is 
consistent  with  the  daytime  visibility  of  2  miles  mea¬ 
sured  earlier  in  the  day.  Those  discrepancies  tha:  occur 
may  have  been  caused  by  the  i  se  of  some  lights  with 
varying  output.  In  actual  practice,  it  would  be  necessary 
to  build  a  library  of  target  lights  with  reasonably  steady 
output.  For  an  initial  sample,  these  results  are  reasonably 
encouraging. 

The  relative  approach  given  by  Eq.  3.19  may  also  be 
illustrated  with  these  data.  Two  examples  are  given  in 
Table  4.  Here  Target  5,  which  yielded  a  somcuhat 
inconsistent  visibility  in  Table  3,  using  the  absolute 
approach,  yields  a  slightly  more  inconsistent  value  using 
the  relative  approach.  Target  6,  which  yielded  a  consis¬ 
tent  result  using  the  absolute  approach  in  Table  3,  also 
yielded  a  consistent  result  using  the  relative  approach  in 
Table  4. 

In  summary,  the  equations  developed  in  Section  3.2 
yielded  reasonable  results  in  this  case  study.  During  the 
day,  the  visibilities  were  around  26  miles  on  16  January 
and  2  miles  on  23  January.  If  one  assumes  that  the 
visibility  (Koschmieder  visibility)  was  still  about  26 
miles  at  night  on  1 6  January,  the  computed  visibility  for 
23  January  at  night  was  about  3  miles.  Given  that  the 
aerosol  may  have  in  fact  changed  with  the  onset  of 
evening,  these  arc  quite  reasonable  results. 

To  automate  this  technique  would  require  several 
steps.  Sensitivity  studies,  to  evaluate  the  relative  merits 
of  the  absolute  and  relative  technique  should  be  done.  In 
addition,  it  would  be  necessary  to  determine  which  light 
sources  are  most  steady,  and  do  a  more  thorough  deter¬ 
mination  of  their  inherent  characteristics.  We  feel  the 
night  visibility  results  are  sufficiently  promising  to  rec¬ 
ommend  continued  development  in  this  area. 

4.  DAYTIME  HSIAVSI  SYSTEM  REFINEMENTS 

A  number  of  improvements  to  the  daytime  capabili¬ 
ties  of  the  HSl/WSI  system  were  evaluated  and/or  in¬ 


stalled.  In  Uiis  next  section,  wc  first  review  the  refine¬ 
ments  which  have  been  fully  programmed.  There  is  a 
discussion  of  an  in-depth  .sensitivity  .study  of  the  H.S1, 
which  revealed  a  number  of  areas  of  potential  system 
improvement.  Finally,  there  is  a  discussion  of  a  recom¬ 
mended  change  to  the  day  visibility  equations  in  use  in  the 
HSI. 

4.1  Changes  to  the  HSI/WSI  Programs 

In  addition  to  the  changes  discussed  in  Section  3.4.  to 
enable  using  the  daytime  program.s  with  the  VS- 100 
imaging  board,  several  other  changes  to  the  daytime 
programs  were  made. 

Pcrltaps  the  most  significant  is  that  the  number  of 
allowable  targets  within  a  scene  was  increased  to  1 6.  This 
is  important  in  increasing  the  accuracy  of  the  determina¬ 
tion,  for  the  following  rca.sons.  As  noted  in  Shield.s.ct  a), 
( 1991 ),  the  error  in  visibility  dclcnminalion  due  to  changes 
in  the  inherent  contrast  of  the  target  with  respect  to  the 
background  can  become  quite  large  when  the  target  range 
becomes  short  with  respect  to  the  visibility,  i.c.  when  V/ 
r  becomes  high.  On  the  other  hand,  the  error  in  vi.sibility 
determination  due  toerrors  in  measured  target  orhorizon 
radiance  become  large  when  the  target  range  approaches 
the  visibility  too  closely,  i.c.  when  V/r  is  loo  close  to  1, 

For  lhc.se  rca.sons,  the  H.Sl  algorithm  has  automatic 
cutoff  values  for  V/r.  If  V/r  lor  a  given  target  is  bctuccn 
1. 15  and  4.0,  the  visibility  determined  using  that  target  is 
given  preference  in  the  sorting  algorithm  which  deter¬ 
mines  the  final  visibility  for  a  given  .sector.  Thus,  on  a 
clear  day,  the  optimum  targets  may  be  quite  distant, 
whereas  on  a  hazy  day,  the  optimum  targets  will  be  closer 
to  the  observer. 

By  allowing  the  u.sc  of  up  to  16  targets  within  a  given 
scene,  it  becomes  possible  to  have  scrs'cral  targets  at  a 
wide  variety  of  ranges.  Thus  it  becomes  more  likely  that 
fora  given  visibility  condition  the  HSIha.s  targets  within 
the  optimal  range,  and  can  thus  yield  a  more  accurate 
value. 


Tabic  4 

Computation  of  23  Jan  Visibility  Using  Reiative  Equation 
Night  Visibility  using  Relative  Approach 


Point  1 

Point  2 

m 

rz 

iNr 

2Nr 

iN'o 

2N„ 

Vis 

(mi) 

Tgl.  2 

2.x 

l.i 

15.6 

27 

!I7 

X9 

6.2 

Effli 

WBsm 

2.8 

1.1 

12.5 

27 

225 

X9 

3.0 

IX 


The  HSI  programs  have  also  been  changed  to  allow 
selection  of  larger  horizon  regions  of  interest  (ROi).  One 
of  the  uncertainties  with  the  HSI  is  the  errorcaused  when 
there  are  clouds  present  in  the  horizon  ROI.  By  selecting 
larger  ROl’s,  the  user  can  evaluate  the  variance  within 
the  horizon  ROI.  We  recommend  using  two  horizon 
ROI’s.and  adding  an  automated  check  which  determines 
if  both  regions  have  similar  averages  and  low  STD’s,  as 
a  means  of  detecting  significant  cloud  contamination. 

Several  additional  program  changes  were  made  for 
user  convenience.  Thc.se  include  saving  a  raw  image 
without  the  ROI’s  overlaid  (in  addition  to  .saving  the 
image  with  ROi’s),  embedding  the  target  information  in 
an  image  header,  saving  the  VVSI  information  todisk,  and 
updating  the  Program  SNGLVIS,  which  allows  later 
evaluation  of  archived  images. 

4.2  Sensitivity  Study  of  Daytime  Visibility 

A  sensitivity  study  was  made,  to  determine  the  most 
significant  sources  of  error  in  the  HSI  system .  This  study. 
Shields,  et  al,  ( 1 99 1 ),  evaluated  the  impact  of  changes  in 
the  inherent  conira.si  of  the  target,  measurement  errors, 
and  changes  in  camera  sensitivity.  The  summary  and 
recommendations  from  this  report  arc  extracted  below. 

4.2.1  Summary  of  the  Sensitivity  Study 

Two  summary  plots,  containing  curves  extracted  from 
the  earlier  sections,  have  been  created:  one  for  Co  =  .8, 
Lq  =  2(X),  and  one  for  Co  =  .5,  Lq  =  2(X).  In  these  plots, 
shown  in  Fig.  15and  16,lhc  fourcurvesshow:  the  impact 
of  a  Co  change  of  .1;  the  impact  of  a  measured  target 
radiance  change  of  4  (on  the  0  to  2.55  scale);  the  impact 
of  a  measured  horizon  radiance  change  of  4;  and  die 
impact  of  system  non-linearity  for  the  system  at  Otis. 

All  of  these  uncertainties  can  cause  a  certain  amount 
of  error.  The  sensitivity  to  Cq  uncertainty  is  very  small 
when  the  target  range  is  close  to  the  visibility  (near  Cr  = 
.05),  and  fairly  large  when  the  target  is  closer.  Unfortu¬ 
nately  however,  the  measurement  uncertainties  cau.se  the 
most  error  when  the  target  range  is  close  to  the  visibility. 
That  is,  when  the  Cq  impacts  arc  least,  the  measurement 
error  impacts  arc  largest. 

There  arc  some  techniques  for  improving  our  Cq 
estimates,  but  in  the  final  analysis  the  Cq  changes  may  be 
mo.st  difficult  to  handle.  The  measurement  uncertainties 
may  in  many  cases  be  mitigated  by  a  combination  of 
improved  measurements  techniques  and  improved  data 
reduction  techniques.  As  the  magnitude  and/or  impact  of 
measurement  uncertainties  arc  improved,  it  should  be 
possible  to  chose  targets  closer  to  the  visual  threshold. 


which  .should  help  mitigate  the  impact  ol  the  C,,  unccr- 
taintic.s. 

Another  way  to  look  at  the  system  is  to  consider  that 
the  ability  to  determine  visibility  from  targets  ranging 
near  the  visibility  depends  on  the  ability  to  accurately 
determine  the  difference  between  signals  which  arc  quite 
close.  This  requires  precision,  stability,  and  accuracy  in 
the  measurements  acquired  by  the  system.  At  the  other 
extreme,  the  ability  to  determine  visibility  from  targets 
which  arc  at  close  range  (and  which  therefore  have 
apparent  contrast  somewhat  close  to  the  inherent  con- 
lra.st)  depends  on  our  ability  to  accurately  characten/e 
the  Co  values  and  their  fluctuations. 

In  tile  short  mn,  there  arc  obvious  ways  in  which 
improvements  can  be  made  in  terms  of  measurements. 
Thc.se  include  keeping  the  horizon  radiance  near  2(K).  and 
measuring  and  applying  the  non-linearity  correction. 
Otherpotential  improvements  such  as  correcting  for  chip 
non-uniformity,  to  improve  measurement  accuracy,  may 
or  may  not  be  warranted.  There  arc  a  variety  of  tests  to 
help  us  determine  questions  such  as  this  one. 

A  significant  improvement  in  the  current  accuracy 
should  be  readily  realized  as  thc.se  changes  arc  enacted. 
In  the  long  run,  as  solid  state  sensors  improve  in  stability 
and  noise  handling,  we  should  expect  significant  im¬ 
provements  in  the  measurement  capabilities.  These  in 
turn  .should  allow  us  to  make  better  use  of  die  measure¬ 
ment  regimes  in  which  the  sensitivity  to  Co  changes 
becomes  .small. 

As  an  outcome  of  this  study,  there  are  a  number  of 
changes  and/or  tests  which  should  be  considered.  These 
fall  roughly  into  two  categories.  The  first  category  is 
changes  having  to  do  with  measurement  accuracy,  in¬ 
cluding  both  improvements  to  the  measurement  accu¬ 
racy,  and  changes  to  mitigate  the  impact  of  measurement 
inaccuracy.  The  .second  category  is  changes  having  to  do 
with  the  handling  of  non-ideal  measurement  conditions, 
such  as  non-ideal  targets  or  non-ideal  horizon  skies. 
These  two  types  of  improvement  categories  will  be 
discussed  below. 

4.2.2  Recommended  Improvements  Relating  to 
Measurement  Accuracy 

The  first  and  most  obvious  change  recore  mended  is  to 
ensure  that  the  horizon  radiance  is  near  2(X)  counts. 
Whereas  this  docs  not  in  itself  increase  the  measurement 
accuracy,  it  significantly  mitigates  the  effccLs  of  mea¬ 
surement  inaccuracy.  With  the  MPL  unit,  the  honzon 
signal  is  cuiremly  near  UK);  we  need  to  verify  that  the 
camera  response  has  not  become  truncated  due  to  intcr- 
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nal  problems,  and  then  adjust  the  auto-iris  to  yield  a 
horizon  brightness  near  200.  With  the  Otis  unit,  we 
similarly  need  to  determine  the  current  signals  which 
occur  for  the  horizon,  and  optimize  them  as  necessary. 

We  need  to  acquire  a  new  linearity  calibration  for  the 
MPL  unit,  to  determine  how  much  change,  if  any,  has 
occurred  since  the  original  calibration.  It  will  be  some¬ 
what  more  difncull  to  acquire  a  linearity  calibration  for 
the  Otis  unit,  but  this  is  certainly  something  to  do  when 
possible.  I  would  propose  that  we  incorporate  applica¬ 
tion  of  the  linearity  correction  into  a  test  version  of  the 
software.  We  should  devise  a  technique  for  testing  the 
sensor  response  in  the  field.  Installation  of  a  light  trap,  to 
enable  testing  of  the  dark  end  of  the  responsivity  curve, 
could  help  us  begin  testing  the  efficacy  of  this  .sort  of  in- 
situ  procedure. 

Next,  there  arc  several  tests  that  involve  investigation 
of  the  magnitude  of  existing  measurement  errors.  As 
discussed  earlier,  both  the  horizon  and  target  radiances 
are  impacted  by  noise,  changes  in  full  dark,  and  chip  non- 
uniformity.  Any  change  in  full  dark  may  be  treated  as 
part  of  the  sensor  responsivity  change  documented  by  the 
linearity  calibration,  and  need  not  be  treated  separately 
for  now. 

The  magnitude  of  system  noise  averaged  over  the 
horizon  and  target  ROI's  may  be  evaluated  by  grabbing 
several  images  of  the  same  scene  in  close  temporal 
succession,  and  comparing  the  resulting  signals  aver¬ 
aged  over  the  ROI’s.  The  errors  due  to  chip  non- 
uniiormity  may  be  evaluated  using  the  linearity  calibra¬ 
tion  data,  if  one  wants  to  know  the  magnitude  of  non¬ 
uniformity  for  a  given  target  location,  one  extracts  the 
signal  for  tliat  target  ROI  and  for  the  horizon  ROI,  but 
using  a  calibration  image  of  a  uniform  source.  Wc  can 
then  decide  if  cither  system  noi.se  or  chip  non-uniformity 
cause  errors  large  enough  to  require  compensation. 

Since  the  impact  of  measurement  error  is  worst  as  V/ 
r  approaches  the  lower  limit  of  1.15,  and  the  impact  of 
non-idea!  Co  is  worst  as  V/r  approaches  4.0,  there  is 
benefit  in  avoiding  these  limits.  Since  the  optimum  range 
is  a  function  of  the  visibility,  it  is  imponant  to  have 
enough  targets  so  that  there  arc  always  several  within  the 
optimum  range,  for  all  possible  visibility  values. 

Within  the  limits  imposed  by  the  availability  of  real 
world  targets,  an  effort  should  be  made  to  select  many 
targets  over  a  range  of  target  distances.  (For  those 
familiar  with  transmissometers,  it  may  be  helpful  to  note 
that  whereas  a  transmissometer  is  reasonably  accurate 
over  a  range  of  visibility  values  dctcimincd  by  its  base 


length,  the  HSl  is  reasonably  accurate  over  a  ranine  ol 
visibility  values  determined  by  the  target  ranges.  By 
using  very  near  targets  under  low  visibility  condiaoas, 
md  far  targets  under  high  visibility  conditions,  wc  are 
creating  an  impact  similar  to  adjusting  the  transmissoni- 
etcr  base  length.) 

Note  also  that  it  is  very  imponant  that  the  range  to 
these  targets  be  determined  accurately .  since  the  error  in 
visibility  due  to  an  error  in  range  is  directly  proportional 
to  the  range  error.  This  may  involve  driving  out  to  the 
sites,  and  visually  identifying  the  targets  being  used, 

4.2.3  Rceommended  Improvcmenis  Relating  to 
Non-idcal  Conditions 

There  are  .several  thing;,  which  can  be  done  to  improve 
the  system  with  regard  to  non-ideal  conditions  First, 
consider  the  impact  of  horizon  sky  problems.  If  the 
measured  horizon  radiance  isnot  equal  to  the  equilibrium 
radiance,  there  is  a  corresponding  error  in  the  visibility. 
Under  clear  sky  conditions,  the  ncar-horizon  radiance  is 
expected  todccrca.se  away  from  the  equilibrium  radiance 
value  as  the  elevation  angle  is  increased.  This  occurs  due 
to  the  decreasing  turbidity  of  the  path  of  .sight. 

It  would  be  instructive  to  extract  the  change  in  horizon 
radiance  over  the  range  of  elevation  angle  in  the  HSi 
ficld-of-vicw.  This  can  be  done  with  existing  HSl 
imagery.  This  would  allow  us  to  determine  the  range  of 
elevation  angles  which  provide  a  surficicntly  accurate 
determinationof  the  equilibrium  radiance  (in  the  absence 
of  clouds). 

Similarly,  wc  should  investigate  the  incidence  of 
clouds  on  the  horizon  which  arc  bnght  enough  to  cause 
error  Wc  can  probably  improve  our  handling  of  this 
po.ssibility  by  using  two  horizon  ROi  ’s,  and  checking  die 
.signal  standard  deviation  in  each,  as  well  as  die  differ¬ 
ence  in  the  average  signals.  If  one  ROI  has  a  high  STD 
and/or  elevated  signal,  use  of  the  other  ROI  might  avoid 
the  cloud.  If  both  ROI's  have  high  STD’s,  there  is 
probably  little  that  can  be  done  currently  other  than  alcn 
the  user,  orpotcntially  not  use  the  visibility  for  that  scene. 
A  test  program  to  investigate  these  possibilities  should  be 
implemented. 

Over  the  longer  run,  the  ca.se  of  cloud  clutter  at  the 
horizon  might  be  addressed  by  the  next  generation  larger 
ficld-of-vicw  WSl.  I f  wc  make  determinations  using  the 
WSl  of  where  the  clouds  arc.  it  should  be  feasible  to 
utilize  horizon  ROI’s  which  arc  in  the  clear  areas  indi¬ 
cated  by  the  WSl.  Similarly  ,  introduction  of  a  rcd/bluc 
filtcrchangcrlothc  HSI  could  allow  determination  of  the 
clear  horizon  regions  for  use  in  visibility  dctcmiination. 
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Finally  comes  the  really  difficult  problem,  changes  is 
Co-  Once  we  have  improved  the  accuracy  of  the  system 
in  the  ways  discussed  above,  we  will  probably  want  to 
tackle  determination  of  Cq.  This  is  done  by  using  the 
median  visibility  determined  from  the  available  targets, 
or  an  independently  determined  visibility,  then  deter¬ 
mining  what  Co  values  for  each  target  will  yield  that 
visibility. 

It  would  be  very  helpful  to  run  curves  similar  to  the 
ones  in  Shields,  el  al,  (1991 ),  which  determine  the  scn.si- 
livity  of  the  Co  determination  to  known  errors,  so  that  we 
can  adjust  our  technique  appropriately.  (For  example, 
we  know  intuitively  that  the  targets  must  be  at  close 
range,  relative  to  the  vi.sibility,  inordcr  for  visibility  to  be 
.sensitive  to  Cq,  which  means  targets  must  be  close  to 
back  out  the  Co  value.  But  does  measurement  error  then 
cause  undue  problems?) 

If  we  arc  successful  with  Cq  extraction,  a  study  of  the 
lime  variation  in  Co  would  be  helpful.  The  equilibrium 
radiance  is  expected  to  change  as  a  function  of  the 
scattering  angle  with  respect  to  the  sun.  This  gives  us  a 
theoretical  change  in  Co  due  to  the  change  in  the  horizon 
brightness.  If  this  is  dominant  over  changes  due  to  target 
brightness,  the  diurnal  changes  in  C©  might  be  reason¬ 
ably  predictable. 

4  3  Recommended  Change  to  the  HSI  Visibility 
Equation 

The  HSI  uses  Eq.  2.2  to  determine  visibility  from  the 
measured  apparent  contrast  of  targets  with  respect  to  the 
horizon.  Al  this  lime,  it  is  recommended  that  the  system 
be  changed  to  use  the  equation 


V  =  r  In  e  /  In 


(4.1) 


The  current  equation  provides  an  estimate  that  is  a 
compromise  between  visibility  and  the  human  estimate 
of  visibility;  this  estimate  is  target  dependent.  Equation 
4.1  would  provide  visibility  determination  consistent 
with  defined  visibility,  and  yield  a  result  which  is  not 
target  dependent  (in  the  absence  of  measurement  error). 


4.3. 1  Derivation  of  The  Koschmieder  Eiquation 

Todcrive  Koschmicdcr’s  Equation,  this  section  stans 
with  the  definitioas  of  radiance  and  contrast,  then  dts- 
cu.sscs  the  definition  of  visibility,  and  show-  how  this 
leads  to  the  desired  equation  for  visibility. 

Radiance 


The  apparent  radiance  of  a  target,  as  measured  irom 
range  r.  is  given  by 


,L,(0.(|))  = 


where 


tl-r(0.<!»  = 

8 

<l> 

tLo  (8, 4))  = 

TrO)  = 

pLr(8,  $)  = 


lL„(8.<p)T,(0)  pLr(O.O)  (4.2) 


apparent  radiance  of  the  target  (rom  a 
range  r  with  look-angle  0,  (J) 
zenith  angle  of  the  path  of  sight 
azimuth  angle  of  the  path  of  sight 
the  inherent  radir  ice  of  the  target,  i  e. 

the  radiance  measured  from  a  range  0 
the  beam  transmittance  of  the  path  of 
sight 

the  path  radiance  of  the  path  of  sight 


These  terms  are  discussed  in  more  detail  in  Duntley, 
(1957).  The  first  term  represents  that  fraction  of  the 
inherent  radiance  which  reaches  the  observer  at  range  r. 
The  second  term,  the  path  radiance,  represents  the  radi¬ 
ance  from  the  surround,  .scattered  into  the  path  of  sight  by 
the  atmosphere  along  the  path  of  sight.  It  is  the  radiance 
which  would  be  observed  if  a  target  were  totally  black. 

Asdiscus.scd  in  detail  in  Duntley,  (1957),  these  terms 
arc  exact  for  monochromatic  radiance,  however  small 
error  is  introduced  by  using  the  phoiopic  pa.ssband  as.so- 
ciated  with  human  vision.  (The  photopic  has  a  passband 
of  100  nm,  centered  at  555  nm.  'fhe  HSI  utilizes  a  filter 
with  a  passband  of  70  nm,  centered  at  550  nnm.)  For 
convenience,  a  short-hand  representation  forEq.  1,  will 
be  used,  as  given  below. 

tW  =  iLoTj-  +  pLr  .  (4.3) 

Contrast 


In  order  to  justify  the  above  statemenLs,  the  following 
sections  will  first  show  a  derivation  of  a  form  of  the 
Koschmiedcrcqualion,  the  equation  as.sociaied  with  clas¬ 
sically  defined  visibility.  Next  will  be  a  derivation  of 
how  visibility  may  be  determined  from  measurements  of 
apparent  contrast  of  targets  o  f  opportunity  using  the  H  SI . 
Finally,  the  difference  between  the  current  and  proposed 
equations  is  discussed. 


A  term  of  importance  in  visibility  theory  is  the  univer¬ 
sal  contrast  of  a  target  with  respect  to  its  background. 
The  contrast  from  range  0  is  defined  as  llic  inherent 
contrast,  given  below. 

.  (4.4) 

b^o 
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Similarly,  the  apparent  contrast,  measured  from  range  r, 
is  given  by 


C_,  =  iij: — ^  .  (4.S) 

bf'r 

To  begin  the  derivation  of  the  Koschmieder  relation, 
one  substitutes  the  expression  for  apparent  radiance  (Eq. 
43)  into  that  for  apparent  contrast  (Eq.  4.5),  then  simpli¬ 
fies  and  regroups  the  expression,  as  shown  below, 

r  -  ~  bt-r  _  +  p^-r|  ~  [b^-oEr  +  pL,] 


(4.6) 


Secondly,  if  the  target  is  precisely  black,  i.c  with 
radiance  equal  to  0,  then  from  Eq.  4,4  it  may  be  seen  that 
Co  becomes  -1.  Since  the  human  visual  threshold  is 
dependent  only  on  the  magnitude  and  not  the  sign  of  the 
contrast,  a  Co  value  of  1  may  be  substituted  into  Eq.  4.8. 
The  equation  thus  becomes 

Cr  =  Tr  .  (4.9) 


The  Resulting  Koschmieder  Equations 

It  is  useful  to  define  the  iraasmissiviiy  T,  which  is  the 
transmittance  over  unit  distance.  (Douglas  and  Booker. 
1977).  Since  the  I'^.nsmittance  TfOver  range  r  is  related 
to  transmissivity  Thy 

Tt  =  (4,10) 

equation  4.9  becomes 


Substituting  the  expression  for  inherent  contrast  (Eq. 
4.4),  this  then  becomes  the  following  equation  for  appar¬ 
ent  contrast  of  the  target. 

C,  .  (4.7) 

bl-f 

This  is  an  expression  for  the  apparent  contrast  of  any 
target  as  viewed  against  any  background. 

Visibility 

The  WMO  defines  visibility  as  follows;  “Meteoro¬ 
logical  visibility  by  day  is  defined  as  the  greatest  distance 
at  which  a  black  object  of  suitable  dimensions,  situated 
near  the  ground,  can  be  seen  and  recognized,  when 
observed  again.st  a  background  of  fog  or  sky,’’  (Quoted 
from  Gordon  1979).  As  will  be  shown,  under  these 
idealized  conditions,  a  human  determination  of  visibility 
approximates  a  mca.sure  of  the  atmospheric  transmit¬ 
tance. 

Idealized  Conditions 

Now  consider  what  happens  if  the  idealized  condi¬ 
tions  implied  by  the  WMOdefiniiion  aremet.  First,  if  the 
path  of  sight  is  horizontal,  and  optically  uniform,  and  ihc 
background  is  a  clear  sky  horizon  radiance,  then  the 
background  radiance  is  the  equilibrium  radiance  (Duntley, 
1957).  Under  these  conditior  s,  as  discussed  in  Duntley 
I957,the  apparent  background  radiance  equals  the  inher¬ 
ent  background  radiance.  Equation  4.7  then  becomes 

Q  =  Co  Tr  (4.8) 

which  is  a  form  of  the  Koschmieder  equation  applying  to 
targets  of  any  inherent  contrast  (Middleton,  p.  70, 1952). 


Cr  =  (4.11) 

To  relate  this  equation  to  the  psychophysically  de¬ 
fined  visibility,  one  may  note  that  the  human  detection  of 
objects  against  the  horizon  is  dependent  on  the  apparent 
contrast  of  the  object.  That  is,  there  is  a  human  contrast 
threshold  e,  which  depends  on  a  number  of  factors  such 
as  targetsize,  as  discussed  in  Gordon  1979.  If  one  defines 
the  visibility  as  the  range  at  which  the  contrast  of  the 
black  target  against  the  horizon  sky  approaches  this 
contrast  threshold  e,  then  one  may  let  r  go  to  V  in  Eq.  4. 1 1 , 
as  the  apparent  contrast  goes  to  threshold.  The  equation 
then  becomes 

e  =  7'^  (4,12) 

which  is  another  form  of  the  Ko.schmieder  equation 
(Douglas  and  Booker,  p.  4-2,  1977). 

Tlii.s  form  of  the  Koschmieder  equation  may  be  con¬ 
sidered  a  definition  of  visibility.  It  corresponds  with  the 
psychophysical  dcfinilion  given  above,  and  is  the  rela¬ 
tion  classically  used  by  sensor  systems  such  as  point 
scatter  meters  and  nephclomelers.  The  WMO  recom¬ 
mends  using  a  threshold  value  e  of  .05  (Douglas  and 
Booker,  1977). 

4.3.2  Determining  the  Vi.sibility  from  Measured 
Apparent  Contrast 

Visibility  as  a  Function  of  Transmissivity 

Using  Equation  4. 12  as  a  definition  of  visibility,  one 
may  convert  the  cxprcs.sion  by  taking  the  log  of  both 
sides,  to  yield  an  expression  for  visibility 

V  =  Ine/ln  T  ,  (4.13) 
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Measured  Transrriissivity 

The  HSl  is  used  to  measure  the  apparent  contrast  of  a 
near-black  target  with  respect  to  the  horizon.  Tnc  target 
does  not  have  inherent  contrast  of  - 1 ,  so  wc  cannot  make 
the  assumptions  made  earlier.  Consider  Eq.  4.7,  repeated 
below. 


or 

( r 

V  r  In  c  /  In  -7^ 

'  .  (4,20) 

This  shows  how  the  classically  defined  visibility  may  be 
determined  from  the  HSI  measurements. 


4.3.3  Relationship  between  the  Proposed  and 
Current  Equation 


The  background  is  still  the  horizon  (since  the  contra.st  of 
tire  target  with  respect  to  the  horizon  is  used).  But  the 
inherent  contrast  is  not  -1.  Therefore  the  equation  dc- 
rcribing  the  apparent  contrast  measured  by  the  HSl  is 

Cr  =  CoTr  =  CoT^  •  (4.15) 


The  current  equation  may  be  derived  by  using  a 
pseudo  visibility,  that  is  by  replacing  Eq.  4.12  with  the 
equation 

^  =  T- 

(4.2!) 


This  apparent  contrast  is  measured  at  some  range  r, 
not  equal  to  visibility.  The  apparent  contrast  in  general 
docs  not  equal  the  contrast  threshold.  However,  by 
knowing  the  range  and  the  inherent  contrast,  this  mea¬ 
surement  of  the  apparent  contrast  may  be  used  as  a 
measure  of  the  transmissivity,  as  follows ,  From  Eq.  4 . 1 5, 
we  have 


C„ 


T’ 


for  the  HSI  measurement,  from  wluch 


In 


=  r  In  r 


(4.16) 


(4.17) 


which  yields  an  expression  for  the  transmissivity. 


In  r  =  -  In 
r 


(4.18) 


This  shows  how  the  transmissivity  of  the  atmosphere 
may  be  determined  from  the  measured  apparent  contrast 
of  a  non-black  target,  with  respect  to  the  horizon,  at  range 
r. 

Resulting  Expression  for  Visibility  from  Measured  Pa¬ 
rameters 


This  is  equivalent  to  using  Eq.  4.8.  rather  than  4,9.  1  his 
defines  the  range  at  which  the  apparent  contra.si  of  the 
non-black  target  reaches  a  human  thrc.shold  contrast. 

Thus  this  pseudo-visibility  V'  is  the  greatest  distance 
at  which  a  NON-black  object  of  suitable  dimensions, 
situated  near  the  ground,  can  be  seen  and  recognized, 
when  observed  against  a  background  of  fog  or  sky.  It  is 
somewhat  less  than  the  visibility,  depending  on  how  non¬ 
black  the  target  is,  For  a  gray  object  ol  0.8  inherent 
contrast,  V'  is  7%  less  than  the  visibility.  (V'  is  not 
equivalent  to  the  visual  range  to  the  target,  unle.ss  the 
target  is  large  and  seen  against  the  horizon  sky.) 

The  proposed  equation  (Eq.  4. 1  or  4.20)  is  not  target 
dependent  (in  the  absence  of  experimental  error),  and  is 
consistent  with  the  definition  of  visibility.  The  current 
equation  (Eq.  2,2)  represents  what  a  human  might  esti¬ 
mate  for  visibility  if  viewing  large  non-black  targets 
against  the  horizon.  The  proposed  equation  represents 
what  a  human  might  estimate  for  visibility  if  viewing 
large  black  targets  against  the  horizon.  Thus,  if  the  shift 
is  made  to  the  proposed  equation,  the  HSl  will  be  more 
self-consistent,  and  more  consistent  with  visibility  as 
defined  by  the  WMO,  but  less  consistent  with  human 
estimates  made  with  non-ideal  targets. 

5.  RELATED  STUDIES 


Now,  this  mca.surcd  iransmissi vi  ty  m  ay  be  substituted 
directly  into  the  defining  equation  for  visibility,  Eq.  4. 1 3, 
to  yield 

V  -  ^ 

^  f  r-  \ 


(4.19) 


Whereas  the  development  work  on  this  contract  was 
directed  primarily  toward  improvement  of  the  visibility 
capability,  related  work  in  the  cloud  determination  capa¬ 
bility,  funded  under  separate  contract,  resulted  in  devel¬ 
opments  which  arc  directly  applicable  to  the  HSl/WSl 
composite  system.  The  first  of  these  related  efforts  was 
development  of  an  improved  cloud  detection  algorithm 
for  the  WSI  (Koehler,  ct  al,  1991),  and  the  second  was 
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development  of  a  night-capable  WSI  system  (Shields,  et 
al.  1991). 

5.1  Directional  Cloud  Decision  Algorithm 

As  discussed  in  Section  2.2,  the  WSI  data  undergo 
several  processing  steps  to  yield  a  cloud  decision  image, 
in  which  opaque  cloud,  thin  cloud,  and  clear  sky  are 
identified.  Several  calibrations  are  applied  in  this  pro¬ 
cess.  and  then  a  red/blue  ratio  image  is  generated. 

This  ratio  image  is  much  more  appropriate  than  the 
original  radiance  image  for  use  in  cloud  detection,  for  a 
number  of  reasons.  It  is  unrealistic  to  use  a  fixed  radiance 
threshold,  because  the  radiance  of  dark  clouds  can  be 
darker  than  the  radiance  of  the  sky.  A  radiance  threshold 
would  also  be  quite  difficult,  because  the  radiance  of  the 
background  sky  varies  considerably  over  the  image,  and 
also  varies  significantly  at  different  times  of  the  day. 

By  using  a  spectral  ratio,  this  first  normalizes  most  of 
the  clear  sky  variance  over  the  image,  as  well  as  most  of 
the  variance  as  a  function  of  time  of  day.  It  also  allows 
a  clean  distinction  between  the  dark  or  gray  clouds  and 
the  dark  blue  sky.  Most  of  the  early  data  processing  of 
WSI  data  used  a  ratio  threshold  to  identify  the  clouds. 

The  opaque  clouds  are  detected  quite  welt  through  use 
of  a  fixed  ratio  threshold.  We  find  that  the  opaque  clouds 
have  a  well  defined  spectral  signature  (as  measured  at  the 
two  wavelengths  used  by  the  WSI).  which  is  reasonably 
independent  of  the  solar  zenith  angle  and  cloud  zenith 
angle.  The  fixed  ratio  threshold  works  well  under  a 
variety  of  conditions,  for  daytime  use  (except  near  dawn 
and  sunset),  and  is  used  in  both  the  original  algorithm  and 
the  new  directional  algoiithm. 

Thin  clouds  are  better  represented  as  a  perturbation 
with  respect  to  the  background  sky  ratio.  Even  though 
the  background  sky  ratio  is  much  less  variable  than  the 
sky  radiance,  there  remains  some  variance  over  the  sky. 
In  particular,  the  ratio  is  higher  near  the  sun  and  near  the 
horizon. 

Several  contrail  cases  were  studied  in  an  effort  to 
determine  the  relationship  between  thin  clouds  and  their 
background  sky  conditions.  We  found  that  the  ratio  of 
thin  cloud  red/bluc  radiance  ratio  to  the  clear  sky  back¬ 
ground  rcdA)luc  radiance  ratio  remains  fairly  constant 
along  sections  of  contrails  exhibiting  uniform  optical 
properties.  I  f  the  clear  sky  background  ratio  were  known, 
a  thin  cloud  discrimination  could  then  be  based  on 
whetherthe  observed  ratioexcceds  the  background  ratio, 
and  the  degree  of  "thinness”  could  be  estimated  by  the 
observed/background  ratio.  The  thin  cloud  decision 


"problem”  then  becomes  one  of  determining  a  reason¬ 
able  estimate  of  the  clear  sky  background  ratio 

The  red/blue  clear  sky  ratio  distribution  is  influenced 
by  many  factors  the  most  important  of  which  are  me  sol  ar 
zenith  angle,  and  the  haze  features  of  the  atmospheric 
boundary  layer.  Figure  1 7  shows  the  variation  of  the  ratio 
disuibution  as  a  functionof  solar  zenith  angle.  Note  how 
the  clear  sky  background  “whitens"  (r/b  ratio  increases) 
as  the  sun  approaches  llte  horizon.  Figure  18  illustrates 
the  boundary  layer  haze  influence,  witli  the  pristine 
February  Arctic  air  mass  yielding  "bluer”  ratios  than  the 
haze  laden  air  mass  from  July. 


Fig.  17  Clear  Sky  Solar  Zenith  Angle  Dependence: 
1  l/Feb/89  •  Columbia,  MO 


Fig.  18  Clear  Sky  Haz.c  L.aycr  Dependence: 
Columbia.  MO,  1989 


Several  approaches  to  modeling  the  background  clear 
sky  ratio  were  evaluated.  It  was  determined  that  due  to 
variations  in  mixing  layer  height  and  prevailing  airma-ss 
type  at  different  sites,  it  was  more  efficacious  to  use 
measurements  of  the  clear  sky  ratio  over  a  period  of  time 
at  a  given  site,  rather  than  model  the  clear  sky  ratio. 

As  the  visibility  changes,  and  the  acro,sol  changes 
from  day  to  day  or  hour  to  hour,  the  magnitude  of  the 
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clear  sky  ratio  changes.  The  modeling  implied  that  these 
changes  arc  proportional  for  the  full  sky  scene.  That  is. 
a  change  from  a  very  clear  day  to  a  hazy  day  results  in  a 
nearly  constant  fractional  change  in  the  r/b  ratio.  The 
WSI  data  were  used  to  verify  that  this  relationship  holds 
reasonably  well  over  a  wide  range  of  conditions.  This 
result  allowed  us  to  use  a  clear  sky  background  which 
depends  on  the  haze  level,  in  the  new  cloud  algorithm .  In 
effect,  this  means  that  a  sky  which  is  white  due  to  haze, 
as  opposed  to  thin  cloud,  is  identified  as  clear. 

In  summary,  the  directional  algorithm  first  identifies 
the  opaque  clouds  through  a  fixed  ratio  threshold.  Then 
the  thin  clouds  are  identified  with  respect  to  a  back* 
ground  clear  sky  r/b  ratio  which  has  been  normalized  for 
haze  level. 

A  sample  cloud  decision  image  resulting  from  the 
directional  algorithm  is  shown  in  Fig.  19.  A  fixed 
threshold  algorithm  would  have  identified  the  down-sun 
portions  of  this  contrail  as  clear.  The  directional  algo¬ 
rithm  is  quite  effective  in  removing  this  directional  bias. 


Fig.  19  Variable  threshold  cloud  decision  image 

The  directional  algorithm  has  not  been  programmed 
into  the  composite  HSI/WSI  system,  due  to  limited 
funding.  However  the  conversion  of  the  system  to 
include  this  update  is  quite  practical. 

5.2  Development  of  the  Day/Night  WSI 

A  major  accomplishment  of  the  Optical  Systems 
Group  has  been  the  development  of  a  WSI  system 
capable  of  acquisition  of  night  imagery. 


The  Day/Night  WSI  takes  advantage  of  many  of  lire 
technologies  used  in  the  Day  WSI.  The  new  sysicn! 
includes  a  larger  fisheye  lens,  with  high  throughput  and 
fuU  upper  hemisphere  coverage  down  to  90  degrees 
zenith  angle.  This  enables  detection  of  the  vertical 
buildupon  the  horizon,  as  well  as  providing  full  coverage 
for  comparison  with  the  meteorological  observer,  li  is 
also  useful  in  characterizing  the  HSI  backgrounds  with 
respect  to  the  presence  of  cloud  clutter. 

The  Day/Night  WSI  uses  a  slow  scan  Charge  Coupled 
Device  (CCD)  sensor  to  yield  the  additional  sensitivity 
required  for  night  operation.  Spectral  filters  are  used,  as 
with  the  Day  WSI,  at  least  down  to  quarter  moon  condi¬ 
tions.  A  fiberoptic  faper  is  used  in  place  of  the  previously 
used  lens  relay  system,  to  resize  and  relocate  the  image 
plane.  The  Day  WSI'sequatorially  driven  .solar  (xcultor 
is  replaced  in  the  Day/Night  system  with  a  zeniUVa/.i- 
muth  dual  drive  occultor.  in  order  to  more  readily  transi  - 
tion  from  sun  to  moon  occultation  and  adapt  to  use  on 
moving  platforms. 

One  of  the  important  design  criteria  is  the  large  range 
of  flux  levels  the  system  must  be  able  to  deal  with.  Figure 
20  shows  the  naturally  occurring  illuminance  levels 
under  a  variety  of  lighting  conditions.  These  data  arc 
from  the  work  of  Brown  1952,  and  are  consistent  with 
irradiance  measurements  acquired  by  our  group  at  the 
Visibility  Lab  over  a  period  of  several  years.  In  Fig,  20, 
the  daytime  illuminance  conditions  the  Day  WSI  has  had 
to  deal  with  are  shown  in  the  top  two  curves  on  the  right 
side  of  the  plot.  These  represent  clear  to  dark  storm 
conditions  for  sun  zenith  angles  0  to  90  degrees.  The 
Day/Night  WSI  can  acquire  imagery  at  450nm  and 
650nm  down  to  quarter  moon  and  urban  starlight  condi¬ 
tions.  For  rural  starlight  conditioas.  it  is  currently  neces¬ 
sary  to  use  open-hole.  TTris  represents  approximately  a  9 
log  range  of  lighting  conditions.  The  sensor  is  designed 
to  obtain  the  necessary  sensitivity  range  by  using  tlic 
approximately  3-4  log  sensitivity  of  the  camera  chip,  2- 
3  logs  range  from  exposure  control,  and  2-5  logs  range 
through  neutral  density  filler  control. 

Two  sample  night  images  are  shown  in  Fig.  2 1 .  Those 
images  were  acquired  in  San  Diego.  Fig.  21a  was 
acquired  under  starlight  (no  moon)  conditions,  with  a  25 
second  exposure.  Flux  levels  were  quite  adequate.  Fig 
21b  was  acquired  under  moonlight.  Only  the  more  dense 
portions  of  the  cloud  layer  were  visually  detected  at  the 
time. 

One  of  the  Day/Night  WSI  units  has  been  fielded  in 
New  Mexico,  and  is  acquiring  very  high  quality  data. 
The  night  cloud  decision  algorithms  arc  in  development 
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Natural  Mlumlnalion  Levels.  These  measurements,  from  Brown 
352  Illustrate  flux  conditions  the  Day/Night  WSl  should  encounter. 


a)  Starlight  (no  moon)  conditions. 


1 1  Sep  91 , 25  seconds  integration 


■  ' 


b)  Moonlight  conditions. 

18  Sep  91, 30  seconds  integration 


Fig.  21  Sample  Night  Sky  Images  from  the  Day/Night  WSI 

at  this  time.  The  Day/Night  WSI  has  many  applications, 
potentially  including  the  automated  weather  system  for 
which  the  composite  HSI/WSI  system  was  developed. 

6.  RECOMMENDATIONS 

The  Composite  HSIAVSI  is  a  very  capable,  unique 
system  applicable  to  both  research  applications  and  auto¬ 


mated  weather  systems.  In  spite  of  its  high  degree  of 
current  capability,  there  remainseveral  areas  of  develop¬ 
ment  which  hold  significant  potential. 

The  first  potential  development  area  to  be  di.scussed 
below  is  conversion  of  the.se  systems  into  tactical  sys¬ 
tems,  or  more  generally  to  small,  relatively  inexpensive 
systems  which  could  be  applied  by  a  broader  user  base. 
The  second  development  area  has  to  do  with  improve¬ 
ment  of  the  quality  and  capabilities  of  the  current  system. 
These  include  further  development  and  testing  of  the 
night  visibility  capability,  and  improvements  in  both  the 
visibility  and  cloud  algorithms.  The  third  development 
area  involves  combining  both  systems  into  one  optical 
system  capable  of  determining  not  only  cloud  cover, 
visibility,  and  radiance  distribution,  but  directional  ter¬ 
rain  reflectance  as  well.  These  three  areas  of  research  and 
development  are  discussed  in  the  follov'ing  sections. 

6.1  Compact  and  Ruggedized  Systems 

In  order  to  make  the  WSI  and  the  HS!  more  available 
to  both  the  research  community  and  the  tactical  military' 
community,  it  is  highly  desirable  to  develop  versions  of 
these  instruments  which  arc  smaller  and  less  costly.  It  is 
important  that  such  systems  require  minimal  interven¬ 
tion  by  personnel  highly  trained  in  the  engineering  of  the 
systems.  The  WSI  and  HSI  units  have  been  operational 
for  several  years;  development  of  compact,  tactical  ver¬ 
sions  of  these  units  is  both  feasible  and  appropriate  at  this 
time. 

6.1.1  A  Compact  WSI  System 

In  order  to  create  a  compact  version  of  the  WSI,  we 
must  consider  the  sensor,  the  controller,  and  the  power 
requirements.  Up  to  this  point,  the  sensor  has  consisted 
of  a  single  solid  stale  camera  (CID  in  the  Day  WSI,  and 
a  slow  scan  CCD  in  the  Day/Nighl  WSI),  with  a  filter 
changer,  aperture  control,  fisheye  and  relay  optics,  and 
occultor. 

It  should  be  feasible  to  miniaturize  this  combination 
primarily  through  use  of  a  color  camera,  to  acquire  the 
necessary  imagery.  In  order  to  evaluate  this  possibility, 
the  Optical  Systems  Group  has  performed  tests  of  a  3- 
chip  color  camera,  acquiring  imagery  of  clouds  and  sky. 
We  determined  that  it  is  possible  to  extract  the  red  and 
blue  components  of  tlw  color  image  separately,  and  ratio 
them  just  as  one  would  the  red  and  blue  images  acquired 
with  the  current  image. 

If  the  Compact  WSI  makes  use  of  such  acolorcamera, 
it  should  not  be  necessary  to  use  a  filter  changer.  Addi¬ 
tionally,  since  the  color  images  would  be  acquired  simul- 


taneously,  rather  than  in  separate  grabs  2  seconds  apart, 
use  of  an  auto  iris  becomes  practical  in  the  compact  WSI. 
The  color  camera  does  not  have  the  same  flexibility  in 
control  of  the  filter  choices,  however  it  has  the  advantage 
of  compact  size.  Also,  due  to  simultaneity  of  the  color 
image  acquisitions,  applications  in  quickly  changing 
environments,  such  as  onboard  an  aircraft,  become  more 
practical. 

There  are  a  variety  of  camera  choices,  depending  on 
the  desired  trade-off  between  resolution  and  cost.  The 
color  camera  used  in  the  above  tests  yielded  full  512  x 
512  images  independently  in  each  color  (RGB),  to  yield 
resolution  equivalent  to  that  of  the  current  WSI.  Less 
costly  cameras  are  available  at  lower  resolution. 

Without  the  need  for  spectral  filters  and  computer 
control  of  the  aperture,  much  of  the  requirement  for 
interface  between  the  computer  and  sensor  is  relieved.  If 
in  addition  the  occultor  is  independently  controlled 
through  a  small  programmable  device  (micro  controller), 
this  eliminates  the  need  for  much  of  the  logic  and  me¬ 
chanical  control  inherent  in  the  current  system. 

It  is  possible  to  further  decrease  the  size  of  the  control 
hardware  by  taking  advantage  of  the  significant  progress 
which  has  been  made  in  recent  years  in  combining  a 
variety  of  functions  on  a  single  computer  board.  With  a 
configuration  much  like  the  current  lap-top  computers, 
possibly  with  a  single  image  board  added,  it  should  be 
possible  to  enable  the  most  basic  capabilities  of  the 
current  WSI  with  a  much  smaller  computer  system. 

Cooling  remains  in  issue  for  the  compact  WSI,  how¬ 
ever  with  a  much  smaller  sensor  package  to  keep  cooled, 
it  should  be  feasible  to  convert  from  the  large  chiller 
system  currently  in  use  to  a  smaller  cooler  such  as  a 
thenno-electric  system.  This  conversion  of  the  cooling 
system  would  also  minimize  the  power  requirements  of 
the  fuU  system. 

The  WSI  system  is  currently  quite  intensive  in  its 
requirement  for  interaction  of  trained  personnel.  This 
interaction  ranges  from  radiometric  calibration,  and  analy¬ 
sis  of  the  calibration  results,  to  evaluation  of  the  clear  sky 
background  used  in  the  directional  algorithm.  Regarding 
the  calibrations,  the  primary  calibrations  involved  are 
image  size  and  flux  control  calibrations,  linearity,  and 
red/blue  relative  radiance  calibrations. 

With  the  compact  senses  package  discussed  above, 
the  first  two  of  these  calibrations  would  no  longer  be 
required.  An  automatic  in-field  calibrator  would  be 
developed  to  test  the  linearity  and  relative  spectral  re¬ 
sponse.  Following  initial  calibration  and  setup,  the 


calibration  could  then  be  automatically  updated  in  the 
field  at  3-raonth  intervals  (or  as  required)  with  the  in¬ 
field  calibrator. 

Regarding  evaluation  of  the  clear  sky  background, 
this  evaluation  is  necessary  only  for  the  thin  cloud 
algorithm.  As  a  first  approach,  the  system  could  be 
trained  to  identify  only  opaque  clouds.  For  thin  clouds, 
if  one  can  accept  lower  resolution  than  that  used  in  the 
current  WSI  system,  there  are  several  possible  algo¬ 
rithms  which  may  prove  to  yield  adequate  thin  cloud 
detection. 

Thus  both  the  WSI  hardware  and  software  have  ma¬ 
tured  to  the  point  that  development  of  a  compact,  more 
transportable  system  is  very  feasible.  The  tactical  mili¬ 
tary  application  is  an  obvious  candidate  for  use  of  such  a 
system.  Likewise,  availability  of  these  units  for  a  variety 
of  research  application  would  be  highly  desirable. 

6.1.2  A  Compact  HSI  System 

The  development  of  a  compact  HSI  unit  would  appear 
to  require  more  technical  compromises  than  required  for 
development  of  a  compact  WSI  unit.  The  scheme  cur¬ 
rently  used  in  the  HSI,  measurement  of  the  apparent 
contrast  of  a  dark  target  with  respect  to  the  sky,  requires 
a  high  resolution  image.  The  HSI  uses  a  very  high 
angular  resolution  rotary  table,  in  order  to  access  the  full 
360  degree  surrounding  horizon  with  the  small  field  of 
view  optics,  acquiring  many  images  at  representative 
azimuthal  points.  The  rotary  table  is  probably  not  prac¬ 
tical  to  compaitmentize  at  this  point. 

A  reasonable  technical  compromise  between  the  need 
to  eliminate  the  rotary  table  and  the  need  for  tight  angular 
control  is  the  use  of  a  hand  held,  narrow  field  of  view 
optical  system  with  crosshairs.  A  user  would  point  the 
crosshairs  at  a  dark  target  of  interest,  acquire  the  image, 
and  then  acquire  the  image  of  the  horizon  sky.  The 
system  could  then  automatically  compute  the  visibility  in 
that  azimuthal  direction.  This  would  require  that  the  user 
input  an  estimate  of  the  range  to  the  target. 

With  this  approach,  the  automation  of  the  HSI  is 
sacrificed  in  return  for  a  convenient,  transportable  unit. 
The  resulting  visibilities  should  be  reasonably  accurate, 
if  an  accurate  range  to  the  target  is  input. 

A  second  approach  to  determination  of  visibility  is  to 
use  the  same  optical  .system  used  for  the  WSI  (either  the 
current  or  the  compact).  The  shape  and  magnitude  of  the 
clear  sky  background  ratio  provides  a  great  deal  of 
information  related  to  the  visibility  and  height  of  the 
inversion  layer.  The  current  cloud  algorithm  makes  use 
of  the  nearly  proportional  relation  between  visibility  and 
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the  magnitude  of  the  clear  sky  background  ratio  at  a  given 
site.  Our  group,  previously  at  the  Visibility  Lab,  has  done 
experimental  work  with  this  inversion  technique  in  the 
past;  sec  Gordon  (1985). 

Thus  a  tactical  visibility  .sensor  seems  rca.sonabIc  to 
evaluate,  however  it  may  not  be  possible  to  retain  the 
automated  capability  in  a  small  tactical  visibility  device 
based  on  the  HSI. 

6.2  Improvements  of  Quality  and  Capability  of 
the  Current  HSI/WSI 

Even  if  development  of  a  compact  WSl  and/or  MSI 
system  comes  to  fruition,  the  current  system  remains  a 
valuable  complement  as  a  high  accuracy  rc.scarch  instru¬ 
ment.  In  this  mode,  the  instrument  may  be  used  for 
making  preci.sion  me  asurcmcntsformodcl  evaluation,  as 
well  as  in  support  of  test  sites.  There  are  several  areas  of 
potential  improvements,  which  will  only  be  summarized 
here,  since  they  have  been  discussed  earlier  in  this  report. 

6.2.1  Integration  of  the  Sensitivity  Study  Results 

As  discussed  in  Section  4.2,  a  sensitivity  study  re¬ 
vealed  several  areas  of  potential  improvement  to  the 
system.  Of  these,  the  most  significant  were  probably  as 
follows.  The  linearity  of  the  camera  should  be  measured 
and  applied  to  the  imagery  on  a  routine  basis.  The 
horizon  region  of  interest  should  be  evaluated  for  cloud 
clutter,  so  that  clear  horizon  is  used.  And  a  better  range 
of  targets  should  be  used,  with  a  sorting  algorithm  that 
gives  more  preference  to  those  targcLs  expected  to  yield 
accurate  results. 

6.2.2  Addition  of  the  Directional  Algorithm 

As  discussed  in  Section  5.1,  a  more  accurate  direc¬ 
tional  cloud  algorithm  has  been  developed  under  sepa¬ 
rate  funding,  which  applies  to  WSI  data.  The  current 
composite  system  contains  the  older  code,  based  on  a 
fixed  ratio  threshold.  This  technique  is  quite  accurate  for 
opaque  clouds  (if  the  proper  calibration  inputs  are  up  to 
date),  however  it  tends  to  under-estimate  the  down-sun 
thin  cloud  amount.  Up  to  this  point,  the  directional  cloud 
algorithm  has  required  some  analyst  intervention  in  the 
selection  of  the  haze  normalization.  As  this  selection 
process  becomes  more  automated,  it  becomes  appropri¬ 
ate  to  convert  the  logic  in  the  composite  system  to  include 
the  directional  cloud  algorithm,  thus  resulting  in  more 
accurate  thin  cloud  assessment. 

6.2.3  Further  Development  of  Night  Visibility 

During  this  year,  theoretical  analysis  was  developed 
which  showed  how  to  determine  visibility  from  measure¬ 


ments  of  lights  of  opponunily  at  night  (Section  3.2). 
Programs  were  developed  to  allow  the  acquisition  of  test 
data  (Section  3.4),  and  the  data  yielded  rea.sonable  results 
(Section  .3.6).  However  an  evaluation  of  the  camera 
{K'rformancc  indicated  that  the  camera  currently  used  in 
the  system  should  tx;  replaced  for  night  applications 
(Section  3.5). 

Wc  recommend  tliat  an  evaluation  ol  alternate  camera 
sen.sorchoiccsshould  be  made,  after  which  furtherdcvcl- 
opment  of  the  automated  night  visibility  capability  can 
proceed.  This  development  should  include  an  error 
analysis  of  the  proposed  theoretical  methods,  program¬ 
ming  of  the  selected  method,  and  evaluation  of  engineer¬ 
ing  considerations  such  as  the  stability  of  available  light 
sources. 

6.3  A  220  Degree  Whole  Sky  Imager 

A  super  wide  angle  version  of  the  WSl  that  presently 
is  under  consideration  at  MPL  would  use  a  220  degree 
field  of  view  super  fisheye  lens  which  is  currently  avail¬ 
able.  Thus  the  full  scene  would  include  the  overhead  sky 
dome,  the  local  horizon,  and  the  surrounding  terrain. 
Output  products  related  to  cloud  cover  would  be  pro¬ 
duced  in  a  manner  similar  to  that  used  in  the  current  WSl 
systems. 

In  addition,  bi-directional  surface  reflectance  data 
may  be  derived  from  the  composite  sun-sky-surfacc 
radiance  measurements.  Ourcxpcricnce  to  date  with  the 
Nikkor  8  mm  fisheye  used  in  the  Day/Night  WSl  (EO 
Camera  6)  illustrates  outstanding  image  quality  and 
angular  resolution  right  up  to  the  outside  edge  of  each 
image.  Assuming  equal  quality  in  the  Nikkor  6mm 
supcrfisheye,  the  detection  and  identification  of  horizon 
targets  and  nearby  surface  features  should  enable  the 
extraction  of  the  desired  surface  radiance.  The  512x512 
chip  would  yield  about  0.4  degrees  per  pixel,  and  0.2 
degrees  per  pixel  may  be  obtained  with  the  use  of  a  1024 
X  1024  chip. 

In  addition  to  cloud  cover,  sky  radiance,  surface 
radiance,  and  surface  reflectance  properties,  this  system 
has  potential  for  determination  of  the  visibility.  Initial 
experimental  evaluation  of  proprietary  techniques  for 
zooming  into  selected  sub.sets  of  the  image  has  yielded 
positive  results.  By  zooming  into  selected  regions  near 
the  horizon,  application  of  algorithms  much  like  those 
used  in  the  HSI  may  become  feasible.  Thus  this  220 
degree  system  has  tremendous  potential  for  a  variety  of 
appi  icalions  requiring  characterization  of  the  atmosphere/ 
surface  radiative  and  meteorological  properties. 
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7.  SUMMARY 

The  Automated  Observing  System  combines  the  ca¬ 
pabilities  of  the  Whole  Sky  Imager  for  cloud  field  assess¬ 
ment,  and  the  Horizon  Scanning  Imager  for  determina¬ 
tion  of  the  sector  visibility.  During  the  past  funding 
interval,  much  of  the  effort  on  this  system  was  directed 
toward  development  of  night  visibility  capability.  A 
theoretical  development  of  the  equations  and  analytic 
approaches  to  night  visibility  has  been  completed.  The 
derived  equations  allow  one  to  rigorously  derive  the 
visibility  from  measurements  acquired  at  night  with  the 
system. 

A  variety  of  hardware  and  software  adaptations  have 
been  completed.  These  include  building  the  electronics 
necessary  to  enable  acquisition  of  imagery  with  multiple 
integration  periods.  This  change  required  extensive 
software  changes,  as  well  as  use  of  a  new  image  board. 
As  a  result  of  the  lack  of  downward  compatibility  in  the 
image  boards,  significant  work  to  adapt  the  Day  WSI  and 
HSI  code  to  the  new  image  board  was  also  required. 

Test  and  calibration  data  indicate  that  the  system  is 
working  reasonably  normally,  although  the  aging  char¬ 
acteristics  of  the  CID  camera  appear  to  make  it  unsuitable 
for  night  visibility.  We  recommend  conversion  of  the 
system  to  a  more  sensitive  CCD.  System  perfonnance 
was  adequate,  however,  to  test  the  night  visibility  scheme. 
The  results  of  the  test  were  very  positive,  yielding  appro¬ 
priate  values  for  visibility  at  night.  Additional  work  in 
the  automation  of  this  technique,  as  well  as  error  analysis, 
is  recommended. 

During  this  period,  an  in-depth  sensitivity  analysis  of 
the  daytime  visibility  system  was  completed.  Several 
areas  for  improvement  of  the  day  visibility  determina¬ 
tions  were  identified. 

Our  recommendations  include  further  development 
of  the  above  issues.  In  addition,  we  feel  that  the  system 
has  evolved  to  the  point  that  development  of  a  tactical 
system  for  cloud  cover  and  possibly  visibility  is  both 
feasible  and  appropriate.  Such  a  system  could  be  por¬ 
table  and  compact,  and  run  on  batteries,  to  provide 
measurements  of  opaque  and  thin  cloud  cover. 
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